CDC Daily by State
davegoblue
2/20/2022
Background
This file is designed to use CDC data to assess coronavirus disease burden by state, including creating and analyzing state-level clusters.
Through March 7, 2021, The COVID Tracking Project collected and integrated data on tests, cases, hospitalizations, deaths, and the like by state and date. The latest code for using this data is available in Coronavirus_Statistics_CTP_v004.Rmd.
The COVID Tracking Project suggest that US federal data sources are now sufficiently robust to be used for analyses that previously relied on COVID Tracking Project. This code is an attempt to update modules in Coronavirus_Statistics_CTP_v004.Rmd to leverage US federal data.
The code in this module builds on code available in _v003, with function and mapping files updated:
- Generic_Added_Utility_Functions_202105_v001.R - generic functions that can be used in other areas
- Coronavirus_CDC_Daily_Functions_v001.R - functions specific to coronavirus daily data
Broadly, the CDC data analyzed by this module includes:
- CDC case and death data by state and date are available for download on the CDC website
- CDC hospital data are available for download on the healthdata.gov website
- CDC vaccines data are also available for download on the CDC website:
Functions and Mapping Files
The tidyverse package is loaded and functions are sourced:
# The tidyverse functions are routinely used without package::function format
library(tidyverse)## -- Attaching packages --------------------------------------- tidyverse 1.3.1 --## v ggplot2 3.3.3 v purrr 0.3.4
## v tibble 3.1.1 v dplyr 1.0.6
## v tidyr 1.1.3 v stringr 1.4.0
## v readr 1.4.0 v forcats 0.5.1## -- Conflicts ------------------------------------------ tidyverse_conflicts() --
## x dplyr::filter() masks stats::filter()
## x dplyr::lag() masks stats::lag()library(geofacet)## Warning: package 'geofacet' was built under R version 4.1.2# Functions are available in source file
source("./Generic_Added_Utility_Functions_202105_v001.R")
source("./Coronavirus_CDC_Daily_Functions_v001.R")A series of mapping files are also available to allow for parameterized processing. Mappings include:
- urlMapper - mapping file for urlType and url location to download data
- renMapper - mapping file for renaming of variables in the raw data file
- selfListMapper - mapping file for transformations by variable type
- fullListMapper - mapping file for transformations across variable types
- lstComboMapper - mapping file for elements to be combined by data type (most common is to combine NYC and NYS data to NY if the file provides them separately)
- uqMapper - mapping file for fields that should combine to be unique keys for processed files
- lstFilterMapper - mapping file for filtering to subset (most common is to keep 50 states and DC)
- vecSelectMapper - mapping file for variables to keep
- checkControlGroupMapper - mapping file for group_by() of control total checks
- checkControlAggMapper - mapping file for numeric variables for control total checks
- checkSimilarityMapper - mapping file for similarity checks to perform
- plotSimilarityMaooer - mapping file for fields where differences in universe should be plotted
- keyAggMapper - mapping file for the aggregate-level control total checks to perform
- perCapMapper - named vector that drives conversion from original field name to per capita field name
- hhsMapper - named vector that drivers numerical variables to keep (and renaming) from HHS capacity data file
These default parameters are maintained in a separate .R file and can be sourced:
source("./Coronavirus_CDC_Daily_Default_Mappings_v002.R")Example Code Processing
The function is run to download and process the latest CDC case, hospitalization, and death data:
readList <- list("cdcDaily"="./RInputFiles/Coronavirus/CDC_dc_downloaded_220220.csv",
"cdcHosp"="./RInputFiles/Coronavirus/CDC_h_downloaded_220220.csv",
"vax"="./RInputFiles/Coronavirus/vaxData_downloaded_220220.csv"
)
compareList <- list("cdcDaily"=readFromRDS("cdc_daily_220206")$dfRaw$cdcDaily,
"cdcHosp"=readFromRDS("cdc_daily_220206")$dfRaw$cdcHosp,
"vax"=readFromRDS("cdc_daily_220206")$dfRaw$vax
)
cdc_daily_220220 <- readRunCDCDaily(thruLabel="Feb 18, 2022",
downloadTo=lapply(readList, FUN=function(x) if(file.exists(x)) NA else x),
readFrom=readList,
compareFile=compareList,
writeLog=NULL,
useClusters=readFromRDS("cdc_daily_210528")$useClusters,
weightedMeanAggs=c("tcpm7", "tdpm7", "cpm7", "dpm7", "hpm7",
"vxcpm7", "vxcgte65pct"
),
skipAssessmentPlots=FALSE,
brewPalette="Paired"
)##
## -- Column specification --------------------------------------------------------
## cols(
## submission_date = col_character(),
## state = col_character(),
## tot_cases = col_double(),
## conf_cases = col_double(),
## prob_cases = col_double(),
## new_case = col_double(),
## pnew_case = col_double(),
## tot_death = col_double(),
## conf_death = col_double(),
## prob_death = col_double(),
## new_death = col_double(),
## pnew_death = col_double(),
## created_at = col_character(),
## consent_cases = col_character(),
## consent_deaths = col_character()
## )##
## *** File has been checked for uniqueness by: state date
##
##
## Checking for similarity of: column names
## In reference but not in current:
## In current but not in reference:
##
## Checking for similarity of: date
## In reference but not in current: 0
## In current but not in reference: 14
##
## Checking for similarity of: state
## In reference but not in current:
## In current but not in reference:
##
##
## ***Differences of at least 5 and at least 5%
##
## date name newValue refValue absDelta pctDelta
## 1 2022-01-30 new_deaths 796 539 257 0.38501873
## 2 2022-01-29 new_deaths 1394 1098 296 0.23756019
## 3 2022-01-23 new_deaths 868 709 159 0.20164870
## 4 2022-01-22 new_deaths 1176 1028 148 0.13430127
## 5 2022-01-16 new_deaths 807 747 60 0.07722008
## 6 2022-01-25 new_deaths 3445 3220 225 0.06751688
## 7 2022-01-24 new_deaths 2679 2505 174 0.06712963
## 8 2022-01-27 new_deaths 2757 2592 165 0.06169377
## 9 2022-01-17 new_deaths 1429 1350 79 0.05685498
## 10 2022-01-26 new_deaths 3023 2858 165 0.05611291
## 11 2022-01-29 new_cases 195076 173891 21185 0.11483412
## 12 2022-01-30 new_cases 138089 124992 13097 0.09956629
## 13 2022-01-31 new_cases 620416 661083 40667 0.06346786
## 14 2022-02-04 new_cases 272825 289747 16922 0.06015941
##
##
## ***Differences of at least 0 and at least 0.1%
##
## state name newValue refValue absDelta pctDelta
## 1 KY tot_deaths 4003629 3992287 11342 0.002836948
## 2 AL tot_deaths 5972978 5963555 9423 0.001578850
## 3 NC tot_deaths 6808527 6798521 10006 0.001470708
## 4 FL tot_cases 1290393798 1286243847 4149951 0.003221214
## 5 MD tot_cases 232491171 231793719 697452 0.003004414
## 6 KY tot_cases 252489934 252077588 412346 0.001634453
## 7 FL new_deaths 68042 66007 2035 0.030362032
## 8 KY new_deaths 13402 13063 339 0.025618742
## 9 AL new_deaths 17741 17371 370 0.021075416
## 10 NC new_deaths 21278 21097 181 0.008542773
## 11 RI new_deaths 3358 3354 4 0.001191895
## 12 MD new_cases 984492 961805 22687 0.023312989
## 13 KY new_cases 1208554 1193647 14907 0.012411118
## 14 TN new_cases 1912511 1926401 13890 0.007236425
## 15 FL new_cases 5648704 5629602 19102 0.003387388
## 16 NC new_cases 2478266 2470242 8024 0.003242998
## 17 SC new_cases 1408611 1405271 3340 0.002373945
## 18 RI new_cases 348326 347901 425 0.001220866
## 19 PW new_cases 2498 2495 3 0.001201682
##
##
##
## Raw file for cdcDaily:
## Rows: 45,540
## Columns: 15
## $ date <date> 2021-12-01, 2020-08-17, 2021-05-31, 2021-07-20, 2020-0~
## $ state <chr> "ND", "MD", "CA", "MD", "VT", "IL", "VT", "MS", "NH", "~
## $ tot_cases <dbl> 163565, 100715, 3685032, 464491, 855, 1130917, 1009, 28~
## $ conf_cases <dbl> 135705, NA, 3685032, NA, NA, 1130917, NA, 176228, NA, 7~
## $ prob_cases <dbl> 27860, NA, 0, NA, NA, 0, NA, 103954, NA, 108997, 0, NA,~
## $ new_cases <dbl> 589, 503, 644, 155, 2, 2304, 10, 1059, 89, 1946, 180, 5~
## $ pnew_case <dbl> 220, 0, 0, 0, 0, 0, 0, 559, 0, 443, 0, 0, 0, 0, NA, 0, ~
## $ tot_deaths <dbl> 1907, 3765, 62011, 9822, 52, 21336, 54, 6730, 86, 12408~
## $ conf_death <dbl> NA, 3616, 62011, 9604, NA, 19306, NA, 4739, NA, 10976, ~
## $ prob_death <dbl> NA, 149, 0, 218, NA, 2030, NA, 1991, NA, 1432, NA, 416,~
## $ new_deaths <dbl> 9, 3, 5, 3, 0, 63, 0, 13, 2, 17, 0, 6, 0, -1, 0, 0, 8, ~
## $ pnew_death <dbl> 0, 0, 0, 1, 0, 16, 0, 7, 0, 2, 0, 0, 0, 0, NA, 0, 0, 4,~
## $ created_at <chr> "12/02/2021 02:35:20 PM", "08/19/2020 12:00:00 AM", "06~
## $ consent_cases <chr> "Agree", "N/A", "Agree", "N/A", "Not agree", "Agree", "~
## $ consent_deaths <chr> "Not agree", "Agree", "Agree", "Agree", "Not agree", "A~##
## -- Column specification --------------------------------------------------------
## cols(
## .default = col_double(),
## state = col_character(),
## date = col_date(format = ""),
## geocoded_state = col_logical()
## )
## i Use `spec()` for the full column specifications.
##
## *** File has been checked for uniqueness by: state date
##
##
## Checking for similarity of: column names
## In reference but not in current:
## In current but not in reference:
##
## Checking for similarity of: date
## In reference but not in current: 0
## In current but not in reference: 15
##
## Checking for similarity of: state
## In reference but not in current:
## In current but not in reference:
##
##
## ***Differences of at least 5 and at least 5%
##
## date name newValue refValue absDelta pctDelta
## 1 2022-02-05 inp 108309 114478 6169 0.05538025
## 2 2022-02-05 hosp_ped 3323 3585 262 0.07585408
## 3 2021-11-24 hosp_ped 1387 1306 81 0.06015596
## 4 2022-02-05 hosp_adult 104794 110893 6099 0.05655417
##
##
## ***Differences of at least 0 and at least 0.1%
##
## state name newValue refValue absDelta pctDelta
## 1 NH hosp_ped 725 811 86 0.111979167
## 2 ME hosp_ped 1373 1431 58 0.041369472
## 3 WV hosp_ped 4435 4554 119 0.026476805
## 4 VT hosp_ped 348 357 9 0.025531915
## 5 AR hosp_ped 10602 10393 209 0.019909502
## 6 KS hosp_ped 3856 3929 73 0.018754014
## 7 SC hosp_ped 7275 7393 118 0.016089446
## 8 VI hosp_ped 81 80 1 0.012422360
## 9 MA hosp_ped 9296 9412 116 0.012401112
## 10 ID hosp_ped 3155 3120 35 0.011155378
## 11 KY hosp_ped 15228 15375 147 0.009606901
## 12 NJ hosp_ped 15981 15838 143 0.008988340
## 13 IN hosp_ped 14697 14787 90 0.006105006
## 14 UT hosp_ped 7026 6998 28 0.003993155
## 15 NV hosp_ped 3856 3871 15 0.003882490
## 16 ND hosp_ped 2898 2909 11 0.003788531
## 17 TN hosp_ped 17497 17563 66 0.003764974
## 18 AL hosp_ped 17263 17319 56 0.003238679
## 19 NC hosp_ped 23574 23649 75 0.003176418
## 20 OR hosp_ped 7333 7356 23 0.003131595
## 21 MO hosp_ped 31461 31363 98 0.003119827
## 22 MS hosp_ped 8953 8926 27 0.003020303
## 23 PA hosp_ped 43632 43509 123 0.002823011
## 24 GA hosp_ped 42185 42079 106 0.002515902
## 25 IA hosp_ped 6153 6168 15 0.002434867
## 26 HI hosp_ped 1909 1913 4 0.002093145
## 27 AZ hosp_ped 22800 22847 47 0.002059281
## 28 NE hosp_ped 6181 6170 11 0.001781232
## 29 WA hosp_ped 10469 10484 15 0.001431776
## 30 CO hosp_ped 17474 17499 25 0.001429674
## 31 WI hosp_ped 8578 8568 10 0.001166453
## 32 IL hosp_ped 35034 35073 39 0.001112585
## 33 OK hosp_ped 20546 20524 22 0.001071342
## 34 RI hosp_ped 2843 2846 3 0.001054667
## 35 PR hosp_ped 16962 16979 17 0.001001738
## 36 AK hosp_ped 1996 1998 2 0.001001502
##
##
##
## Raw file for cdcHosp:
## Rows: 38,675
## Columns: 117
## $ state <chr> ~
## $ date <date> ~
## $ critical_staffing_shortage_today_yes <dbl> ~
## $ critical_staffing_shortage_today_no <dbl> ~
## $ critical_staffing_shortage_today_not_reported <dbl> ~
## $ critical_staffing_shortage_anticipated_within_week_yes <dbl> ~
## $ critical_staffing_shortage_anticipated_within_week_no <dbl> ~
## $ critical_staffing_shortage_anticipated_within_week_not_reported <dbl> ~
## $ hospital_onset_covid <dbl> ~
## $ hospital_onset_covid_coverage <dbl> ~
## $ inpatient_beds <dbl> ~
## $ inpatient_beds_coverage <dbl> ~
## $ inpatient_beds_used <dbl> ~
## $ inpatient_beds_used_coverage <dbl> ~
## $ inp <dbl> ~
## $ inpatient_beds_used_covid_coverage <dbl> ~
## $ previous_day_admission_adult_covid_confirmed <dbl> ~
## $ previous_day_admission_adult_covid_confirmed_coverage <dbl> ~
## $ previous_day_admission_adult_covid_suspected <dbl> ~
## $ previous_day_admission_adult_covid_suspected_coverage <dbl> ~
## $ previous_day_admission_pediatric_covid_confirmed <dbl> ~
## $ previous_day_admission_pediatric_covid_confirmed_coverage <dbl> ~
## $ previous_day_admission_pediatric_covid_suspected <dbl> ~
## $ previous_day_admission_pediatric_covid_suspected_coverage <dbl> ~
## $ staffed_adult_icu_bed_occupancy <dbl> ~
## $ staffed_adult_icu_bed_occupancy_coverage <dbl> ~
## $ staffed_icu_adult_patients_confirmed_and_suspected_covid <dbl> ~
## $ staffed_icu_adult_patients_confirmed_and_suspected_covid_coverage <dbl> ~
## $ staffed_icu_adult_patients_confirmed_covid <dbl> ~
## $ staffed_icu_adult_patients_confirmed_covid_coverage <dbl> ~
## $ hosp_adult <dbl> ~
## $ total_adult_patients_hospitalized_confirmed_and_suspected_covid_coverage <dbl> ~
## $ total_adult_patients_hospitalized_confirmed_covid <dbl> ~
## $ total_adult_patients_hospitalized_confirmed_covid_coverage <dbl> ~
## $ hosp_ped <dbl> ~
## $ total_pediatric_patients_hospitalized_confirmed_and_suspected_covid_coverage <dbl> ~
## $ total_pediatric_patients_hospitalized_confirmed_covid <dbl> ~
## $ total_pediatric_patients_hospitalized_confirmed_covid_coverage <dbl> ~
## $ total_staffed_adult_icu_beds <dbl> ~
## $ total_staffed_adult_icu_beds_coverage <dbl> ~
## $ inpatient_beds_utilization <dbl> ~
## $ inpatient_beds_utilization_coverage <dbl> ~
## $ inpatient_beds_utilization_numerator <dbl> ~
## $ inpatient_beds_utilization_denominator <dbl> ~
## $ percent_of_inpatients_with_covid <dbl> ~
## $ percent_of_inpatients_with_covid_coverage <dbl> ~
## $ percent_of_inpatients_with_covid_numerator <dbl> ~
## $ percent_of_inpatients_with_covid_denominator <dbl> ~
## $ inpatient_bed_covid_utilization <dbl> ~
## $ inpatient_bed_covid_utilization_coverage <dbl> ~
## $ inpatient_bed_covid_utilization_numerator <dbl> ~
## $ inpatient_bed_covid_utilization_denominator <dbl> ~
## $ adult_icu_bed_covid_utilization <dbl> ~
## $ adult_icu_bed_covid_utilization_coverage <dbl> ~
## $ adult_icu_bed_covid_utilization_numerator <dbl> ~
## $ adult_icu_bed_covid_utilization_denominator <dbl> ~
## $ adult_icu_bed_utilization <dbl> ~
## $ adult_icu_bed_utilization_coverage <dbl> ~
## $ adult_icu_bed_utilization_numerator <dbl> ~
## $ adult_icu_bed_utilization_denominator <dbl> ~
## $ geocoded_state <lgl> ~
## $ `previous_day_admission_adult_covid_confirmed_18-19` <dbl> ~
## $ `previous_day_admission_adult_covid_confirmed_18-19_coverage` <dbl> ~
## $ `previous_day_admission_adult_covid_confirmed_20-29` <dbl> ~
## $ `previous_day_admission_adult_covid_confirmed_20-29_coverage` <dbl> ~
## $ `previous_day_admission_adult_covid_confirmed_30-39` <dbl> ~
## $ `previous_day_admission_adult_covid_confirmed_30-39_coverage` <dbl> ~
## $ `previous_day_admission_adult_covid_confirmed_40-49` <dbl> ~
## $ `previous_day_admission_adult_covid_confirmed_40-49_coverage` <dbl> ~
## $ `previous_day_admission_adult_covid_confirmed_50-59` <dbl> ~
## $ `previous_day_admission_adult_covid_confirmed_50-59_coverage` <dbl> ~
## $ `previous_day_admission_adult_covid_confirmed_60-69` <dbl> ~
## $ `previous_day_admission_adult_covid_confirmed_60-69_coverage` <dbl> ~
## $ `previous_day_admission_adult_covid_confirmed_70-79` <dbl> ~
## $ `previous_day_admission_adult_covid_confirmed_70-79_coverage` <dbl> ~
## $ `previous_day_admission_adult_covid_confirmed_80+` <dbl> ~
## $ `previous_day_admission_adult_covid_confirmed_80+_coverage` <dbl> ~
## $ previous_day_admission_adult_covid_confirmed_unknown <dbl> ~
## $ previous_day_admission_adult_covid_confirmed_unknown_coverage <dbl> ~
## $ `previous_day_admission_adult_covid_suspected_18-19` <dbl> ~
## $ `previous_day_admission_adult_covid_suspected_18-19_coverage` <dbl> ~
## $ `previous_day_admission_adult_covid_suspected_20-29` <dbl> ~
## $ `previous_day_admission_adult_covid_suspected_20-29_coverage` <dbl> ~
## $ `previous_day_admission_adult_covid_suspected_30-39` <dbl> ~
## $ `previous_day_admission_adult_covid_suspected_30-39_coverage` <dbl> ~
## $ `previous_day_admission_adult_covid_suspected_40-49` <dbl> ~
## $ `previous_day_admission_adult_covid_suspected_40-49_coverage` <dbl> ~
## $ `previous_day_admission_adult_covid_suspected_50-59` <dbl> ~
## $ `previous_day_admission_adult_covid_suspected_50-59_coverage` <dbl> ~
## $ `previous_day_admission_adult_covid_suspected_60-69` <dbl> ~
## $ `previous_day_admission_adult_covid_suspected_60-69_coverage` <dbl> ~
## $ `previous_day_admission_adult_covid_suspected_70-79` <dbl> ~
## $ `previous_day_admission_adult_covid_suspected_70-79_coverage` <dbl> ~
## $ `previous_day_admission_adult_covid_suspected_80+` <dbl> ~
## $ `previous_day_admission_adult_covid_suspected_80+_coverage` <dbl> ~
## $ previous_day_admission_adult_covid_suspected_unknown <dbl> ~
## $ previous_day_admission_adult_covid_suspected_unknown_coverage <dbl> ~
## $ deaths_covid <dbl> ~
## $ deaths_covid_coverage <dbl> ~
## $ on_hand_supply_therapeutic_a_casirivimab_imdevimab_courses <dbl> ~
## $ on_hand_supply_therapeutic_b_bamlanivimab_courses <dbl> ~
## $ on_hand_supply_therapeutic_c_bamlanivimab_etesevimab_courses <dbl> ~
## $ previous_week_therapeutic_a_casirivimab_imdevimab_courses_used <dbl> ~
## $ previous_week_therapeutic_b_bamlanivimab_courses_used <dbl> ~
## $ previous_week_therapeutic_c_bamlanivimab_etesevimab_courses_used <dbl> ~
## $ icu_patients_confirmed_influenza <dbl> ~
## $ icu_patients_confirmed_influenza_coverage <dbl> ~
## $ previous_day_admission_influenza_confirmed <dbl> ~
## $ previous_day_admission_influenza_confirmed_coverage <dbl> ~
## $ previous_day_deaths_covid_and_influenza <dbl> ~
## $ previous_day_deaths_covid_and_influenza_coverage <dbl> ~
## $ previous_day_deaths_influenza <dbl> ~
## $ previous_day_deaths_influenza_coverage <dbl> ~
## $ total_patients_hospitalized_confirmed_influenza <dbl> ~
## $ total_patients_hospitalized_confirmed_influenza_and_covid <dbl> ~
## $ total_patients_hospitalized_confirmed_influenza_and_covid_coverage <dbl> ~
## $ total_patients_hospitalized_confirmed_influenza_coverage <dbl> ~##
## -- Column specification --------------------------------------------------------
## cols(
## .default = col_double(),
## Date = col_character(),
## Location = col_character()
## )
## i Use `spec()` for the full column specifications.
##
## *** File has been checked for uniqueness by: state date
##
##
## Checking for similarity of: column names
## In reference but not in current:
## In current but not in reference:
##
## Checking for similarity of: date
## In reference but not in current: 0
## In current but not in reference: 14
##
## Checking for similarity of: state
## In reference but not in current:
## In current but not in reference:
##
##
## ***Differences of at least 1 and at least 1%
##
## [1] date name newValue refValue absDelta pctDelta
## <0 rows> (or 0-length row.names)##
##
## ***Differences of at least 0 and at least 0.1%
##
## [1] state name newValue refValue absDelta pctDelta
## <0 rows> (or 0-length row.names)
##
##
##
## Raw file for vax:
## Rows: 27,992
## Columns: 82
## $ date <date> 2022-02-19, 2022-02-19, 2022-0~
## $ MMWR_week <dbl> 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7~
## $ state <chr> "NC", "TN", "MN", "MI", "SD", "~
## $ Distributed <dbl> 20744900, 12186030, 11914970, 1~
## $ Distributed_Janssen <dbl> 916100, 503900, 500200, 926300,~
## $ Distributed_Moderna <dbl> 7813760, 4644240, 4216760, 7835~
## $ Distributed_Pfizer <dbl> 12015040, 7037890, 7198010, 111~
## $ Distributed_Unk_Manuf <dbl> 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0~
## $ Dist_Per_100K <dbl> 197795, 178441, 211272, 199329,~
## $ Distributed_Per_100k_12Plus <dbl> 230870, 208823, 249367, 231608,~
## $ Distributed_Per_100k_18Plus <dbl> 253377, 229098, 274762, 253818,~
## $ Distributed_Per_100k_65Plus <dbl> 1184680, 1065780, 1294570, 1127~
## $ vxa <dbl> 16040239, 9551129, 9853584, 150~
## $ Administered_12Plus <dbl> 15576577, 9369683, 9460637, 146~
## $ Administered_18Plus <dbl> 14630091, 8914389, 8826086, 138~
## $ Administered_65Plus <dbl> 4239236, 2778123, 2487485, 4293~
## $ Administered_Janssen <dbl> 508845, 259901, 353693, 459665,~
## $ Administered_Moderna <dbl> 5969173, 3654211, 3581985, 5900~
## $ Administered_Pfizer <dbl> 9561290, 5583600, 5913885, 8724~
## $ Administered_Unk_Manuf <dbl> 931, 53417, 4021, 2106, 133, 32~
## $ Admin_Per_100k <dbl> 152938, 139858, 174720, 151062,~
## $ Admin_Per_100k_12Plus <dbl> 173352, 160562, 198000, 170666,~
## $ Admin_Per_100k_18Plus <dbl> 178691, 167591, 203531, 176626,~
## $ Admin_Per_100k_65Plus <dbl> 242091, 242972, 270267, 243193,~
## $ Recip_Administered <dbl> 15939232, 9383280, 9868373, 153~
## $ Administered_Dose1_Recip <dbl> 8596653, 4180275, 4183752, 6576~
## $ Administered_Dose1_Pop_Pct <dbl> 82.0, 61.2, 74.2, 65.9, 74.7, 0~
## $ Administered_Dose1_Recip_12Plus <dbl> 8331132, 4079234, 3968078, 6351~
## $ Administered_Dose1_Recip_12PlusPop_Pct <dbl> 92.7, 69.9, 83.0, 73.9, 86.3, 0~
## $ Administered_Dose1_Recip_18Plus <dbl> 7823417, 3850407, 3680383, 5966~
## $ Administered_Dose1_Recip_18PlusPop_Pct <dbl> 95.0, 72.4, 84.9, 76.1, 88.9, 0~
## $ Administered_Dose1_Recip_65Plus <dbl> 2154949, 1047531, 937204, 16884~
## $ Administered_Dose1_Recip_65PlusPop_Pct <dbl> 95.0, 91.6, 95.0, 95.0, 95.0, 0~
## $ vxc <dbl> 6201249, 3646584, 3830382, 5889~
## $ vxcpoppct <dbl> 59.1, 53.4, 67.9, 59.0, 59.7, 0~
## $ Series_Complete_12Plus <dbl> 6011177, 3568185, 3651613, 5701~
## $ Series_Complete_12PlusPop_Pct <dbl> 66.9, 61.1, 76.4, 66.3, 69.0, 0~
## $ vxcgte18 <dbl> 5622542, 3375377, 3382210, 5355~
## $ vxcgte18pct <dbl> 68.7, 63.5, 78.0, 68.3, 71.3, 0~
## $ vxcgte65 <dbl> 1498685, 956337, 876804, 153840~
## $ vxcgte65pct <dbl> 85.6, 83.6, 95.0, 87.1, 92.5, 0~
## $ Series_Complete_Janssen <dbl> 477185, 232189, 326034, 416641,~
## $ Series_Complete_Moderna <dbl> 2152155, 1299423, 1287594, 2139~
## $ Series_Complete_Pfizer <dbl> 3571765, 2103546, 2215307, 3333~
## $ Series_Complete_Unk_Manuf <dbl> 144, 11426, 1447, 1082, 34, 0, ~
## $ Series_Complete_Janssen_12Plus <dbl> 477158, 232135, 326016, 416612,~
## $ Series_Complete_Moderna_12Plus <dbl> 2152040, 1299371, 1287540, 2138~
## $ Series_Complete_Pfizer_12Plus <dbl> 3381836, 2025317, 2036626, 3144~
## $ Series_Complete_Unk_Manuf_12Plus <dbl> 143, 11362, 1431, 1073, 34, 0, ~
## $ Series_Complete_Janssen_18Plus <dbl> 475728, 231891, 325496, 416312,~
## $ Series_Complete_Moderna_18Plus <dbl> 2149019, 1298802, 1285260, 2138~
## $ Series_Complete_Pfizer_18Plus <dbl> 2997656, 1833427, 1770066, 2799~
## $ Series_Complete_Unk_Manuf_18Plus <dbl> 139, 11257, 1388, 986, 34, 0, 5~
## $ Series_Complete_Janssen_65Plus <dbl> 54321, 35691, 50477, 70861, 498~
## $ Series_Complete_Moderna_65Plus <dbl> 720300, 474871, 369087, 768663,~
## $ Series_Complete_Pfizer_65Plus <dbl> 723999, 439839, 456889, 698286,~
## $ Series_Complete_Unk_Manuf_65Plus <dbl> 65, 5936, 351, 592, 21, 0, 2511~
## $ Additional_Doses <dbl> 1544360, 1529958, 2125396, 2985~
## $ Additional_Doses_Vax_Pct <dbl> 24.9, 42.0, 55.5, 50.7, 39.8, 2~
## $ Additional_Doses_12Plus <dbl> 1544252, 1529687, 2125156, 2985~
## $ Additional_Doses_12Plus_Vax_Pct <dbl> 25.7, 42.9, 58.2, 52.4, 41.2, 2~
## $ Additional_Doses_18Plus <dbl> 1500845, 1502838, 2049913, 2906~
## $ Additional_Doses_18Plus_Vax_Pct <dbl> 26.7, 44.5, 60.6, 54.3, 43.1, 2~
## $ Additional_Doses_50Plus <dbl> 1017165, 1059552, 1281534, 1976~
## $ Additional_Doses_50Plus_Vax_Pct <dbl> 33.8, 56.3, 72.6, 64.9, 54.8, 4~
## $ Additional_Doses_65Plus <dbl> 578981, 632802, 708477, 1135879~
## $ Additional_Doses_65Plus_Vax_Pct <dbl> 38.6, 66.2, 80.8, 73.8, 62.9, 5~
## $ Additional_Doses_Moderna <dbl> 680325, 649042, 857058, 1316220~
## $ Additional_Doses_Pfizer <dbl> 836934, 856740, 1239887, 162309~
## $ Additional_Doses_Janssen <dbl> 27082, 20983, 28141, 46218, 254~
## $ Additional_Doses_Unk_Manuf <dbl> 19, 3193, 310, 106, 9, 22, 648,~
## $ Administered_Dose1_Recip_5Plus <dbl> 8594663, 4179589, 4181728, 6576~
## $ Administered_Dose1_Recip_5PlusPop_Pct <dbl> 87.0, 65.1, 79.1, 69.8, 80.3, 0~
## $ Series_Complete_5Plus <dbl> 6200658, 3646444, 3829701, 5889~
## $ Series_Complete_5PlusPop_Pct <dbl> 62.8, 56.8, 72.4, 62.5, 64.1, 0~
## $ Administered_5Plus <dbl> 16037693, 9550281, 9850893, 150~
## $ Admin_Per_100k_5Plus <dbl> 162353, 148745, 186287, 160139,~
## $ Distributed_Per_100k_5Plus <dbl> 210004, 189797, 225320, 211315,~
## $ Series_Complete_Moderna_5Plus <dbl> 2152112, 1299406, 1287586, 2138~
## $ Series_Complete_Pfizer_5Plus <dbl> 3571235, 2103464, 2214653, 3333~
## $ Series_Complete_Janssen_5Plus <dbl> 477168, 232150, 326019, 416627,~
## $ Series_Complete_Unk_Manuf_5Plus <dbl> 143, 11424, 1443, 1081, 34, 0, ~
##
## Column sums before and after applying filtering rules:
## # A tibble: 3 x 6
## isType tot_cases tot_deaths new_cases new_deaths n
## <chr> <dbl> <dbl> <dbl> <dbl> <dbl>
## 1 before 1.84e+10 3.17e+8 7.79e+7 910251 44781
## 2 after 1.83e+10 3.15e+8 7.73e+7 905741 38709
## 3 pctchg 4.83e- 3 4.32e-3 6.93e-3 0.00495 0.136
##
##
## Processed for cdcDaily:
## Rows: 38,709
## Columns: 6
## $ date <date> 2021-12-01, 2020-08-17, 2021-05-31, 2021-07-20, 2020-05-13~
## $ state <chr> "ND", "MD", "CA", "MD", "VT", "IL", "VT", "MS", "NH", "NC",~
## $ tot_cases <dbl> 163565, 100715, 3685032, 464491, 855, 1130917, 1009, 280182~
## $ tot_deaths <dbl> 1907, 3765, 62011, 9822, 52, 21336, 54, 6730, 86, 12408, 55~
## $ new_cases <dbl> 589, 503, 644, 155, 2, 2304, 10, 1059, 89, 1946, 180, 537, ~
## $ new_deaths <dbl> 9, 3, 5, 3, 0, 63, 0, 13, 2, 17, 0, 6, 0, -1, 0, 0, 8, 11, ~
##
## Column sums before and after applying filtering rules:
## # A tibble: 3 x 5
## isType inp hosp_adult hosp_ped n
## <chr> <dbl> <dbl> <dbl> <dbl>
## 1 before 4.51e+7 3.88e+7 945228 38675
## 2 after 4.49e+7 3.86e+7 927959 37083
## 3 pctchg 4.84e-3 4.63e-3 0.0183 0.0412
##
##
## Processed for cdcHosp:
## Rows: 37,083
## Columns: 5
## $ date <date> 2020-10-14, 2020-10-14, 2020-10-11, 2020-10-10, 2020-10-09~
## $ state <chr> "HI", "NE", "IA", "NH", "HI", "DC", "KS", "NM", "ME", "NE",~
## $ inp <dbl> 111, 376, 497, 45, 110, 166, 474, 189, 23, 316, 546, 3246, ~
## $ hosp_adult <dbl> 111, 367, 487, 44, 108, 149, 454, 186, 23, 315, 534, 3104, ~
## $ hosp_ped <dbl> 0, 9, 10, 1, 2, 17, 5, 3, 0, 6, 12, 55, 8, 0, 1, 8, 2, 8, 6~
##
## Column sums before and after applying filtering rules:
## # A tibble: 3 x 9
## isType vxa vxc vxcpoppct vxcgte65 vxcgte65pct vxcgte18 vxcgte18pct
## <chr> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl>
## 1 before 2.66e+11 1.13e+11 1003870. 3.03e+10 1559557. 1.06e+11 1202494.
## 2 after 1.28e+11 5.46e+10 843159. 1.46e+10 1396120. 5.14e+10 1020516.
## 3 pctchg 5.20e- 1 5.16e- 1 0.160 5.16e- 1 0.105 5.17e- 1 0.151
## # ... with 1 more variable: n <dbl>
##
##
## Processed for vax:
## Rows: 22,083
## Columns: 9
## $ date <date> 2022-02-19, 2022-02-19, 2022-02-19, 2022-02-19, 2022-02-1~
## $ state <chr> "NC", "TN", "MN", "MI", "SD", "OH", "MT", "WV", "VA", "IA"~
## $ vxa <dbl> 16040239, 9551129, 9853584, 15086338, 1349798, 17152418, 1~
## $ vxc <dbl> 6201249, 3646584, 3830382, 5889772, 527824, 6712161, 59625~
## $ vxcpoppct <dbl> 59.1, 53.4, 67.9, 59.0, 59.7, 57.4, 55.8, 56.6, 71.7, 60.9~
## $ vxcgte65 <dbl> 1498685, 956337, 876804, 1538402, 140420, 1779459, 175563,~
## $ vxcgte65pct <dbl> 85.6, 83.6, 95.0, 87.1, 92.5, 87.0, 85.0, 83.5, 91.2, 92.0~
## $ vxcgte18 <dbl> 5622542, 3375377, 3382210, 5355218, 476217, 6118597, 54652~
## $ vxcgte18pct <dbl> 68.7, 63.5, 78.0, 68.3, 71.3, 67.2, 65.0, 65.8, 81.0, 71.7~
##
## Integrated per capita data file:
## Rows: 38,973
## Columns: 34
## $ date <date> 2020-01-01, 2020-01-01, 2020-01-01, 2020-01-01, 2020-01-0~
## $ state <chr> "AL", "HI", "IN", "LA", "MN", "MT", "NC", "TX", "AL", "HI"~
## $ tot_cases <dbl> NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA~
## $ tot_deaths <dbl> NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA~
## $ new_cases <dbl> NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA~
## $ new_deaths <dbl> NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA~
## $ inp <dbl> NA, 0, 0, NA, 0, 0, 0, 0, NA, 0, 0, NA, 0, 0, 0, 1877, 0, ~
## $ hosp_adult <dbl> NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA~
## $ hosp_ped <dbl> NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA~
## $ vxa <dbl> NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA~
## $ vxc <dbl> NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA~
## $ vxcpoppct <dbl> NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA~
## $ vxcgte65 <dbl> NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA~
## $ vxcgte65pct <dbl> NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA~
## $ vxcgte18 <dbl> NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA~
## $ vxcgte18pct <dbl> NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA~
## $ tcpm <dbl> NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA~
## $ tdpm <dbl> NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA~
## $ cpm <dbl> NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA~
## $ dpm <dbl> NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA~
## $ hpm <dbl> NA, 0.0000, 0.0000, NA, 0.0000, 0.0000, 0.0000, 0.0000, NA~
## $ ahpm <dbl> NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA~
## $ phpm <dbl> NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA~
## $ vxapm <dbl> NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA~
## $ vxcpm <dbl> NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA~
## $ tcpm7 <dbl> NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA~
## $ tdpm7 <dbl> NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA~
## $ cpm7 <dbl> NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA~
## $ dpm7 <dbl> NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA~
## $ hpm7 <dbl> NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA~
## $ ahpm7 <dbl> NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA~
## $ phpm7 <dbl> NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA~
## $ vxapm7 <dbl> NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA~
## $ vxcpm7 <dbl> NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA~## Warning in showSRID(uprojargs, format = "PROJ", multiline = "NO", prefer_proj =
## prefer_proj): Discarded datum unknown in CRS definition
saveToRDS(cdc_daily_220220, ovrWriteError=FALSE)The latest hospital data are downloaded:
# Run for latest data, save as RDS
indivHosp_20220221 <- downloadReadHospitalData(loc="./RInputFiles/Coronavirus/HHS_Hospital_20220221.csv")##
## -- Column specification --------------------------------------------------------
## cols(
## .default = col_double(),
## hospital_pk = col_character(),
## collection_week = col_date(format = ""),
## state = col_character(),
## ccn = col_character(),
## hospital_name = col_character(),
## address = col_character(),
## city = col_character(),
## zip = col_character(),
## hospital_subtype = col_character(),
## fips_code = col_character(),
## is_metro_micro = col_logical(),
## geocoded_hospital_address = col_character(),
## hhs_ids = col_character(),
## is_corrected = col_logical()
## )
## i Use `spec()` for the full column specifications.## Rows: 399,863
## Columns: 109
## $ hospital_pk <chr> ~
## $ collection_week <date> ~
## $ state <chr> ~
## $ ccn <chr> ~
## $ hospital_name <chr> ~
## $ address <chr> ~
## $ city <chr> ~
## $ zip <chr> ~
## $ hospital_subtype <chr> ~
## $ fips_code <chr> ~
## $ is_metro_micro <lgl> ~
## $ total_beds_7_day_avg <dbl> ~
## $ all_adult_hospital_beds_7_day_avg <dbl> ~
## $ all_adult_hospital_inpatient_beds_7_day_avg <dbl> ~
## $ inpatient_beds_used_7_day_avg <dbl> ~
## $ all_adult_hospital_inpatient_bed_occupied_7_day_avg <dbl> ~
## $ inpatient_beds_used_covid_7_day_avg <dbl> ~
## $ total_adult_patients_hospitalized_confirmed_and_suspected_covid_7_day_avg <dbl> ~
## $ total_adult_patients_hospitalized_confirmed_covid_7_day_avg <dbl> ~
## $ total_pediatric_patients_hospitalized_confirmed_and_suspected_covid_7_day_avg <dbl> ~
## $ total_pediatric_patients_hospitalized_confirmed_covid_7_day_avg <dbl> ~
## $ inpatient_beds_7_day_avg <dbl> ~
## $ total_icu_beds_7_day_avg <dbl> ~
## $ total_staffed_adult_icu_beds_7_day_avg <dbl> ~
## $ icu_beds_used_7_day_avg <dbl> ~
## $ staffed_adult_icu_bed_occupancy_7_day_avg <dbl> ~
## $ staffed_icu_adult_patients_confirmed_and_suspected_covid_7_day_avg <dbl> ~
## $ staffed_icu_adult_patients_confirmed_covid_7_day_avg <dbl> ~
## $ total_patients_hospitalized_confirmed_influenza_7_day_avg <dbl> ~
## $ icu_patients_confirmed_influenza_7_day_avg <dbl> ~
## $ total_patients_hospitalized_confirmed_influenza_and_covid_7_day_avg <dbl> ~
## $ total_beds_7_day_sum <dbl> ~
## $ all_adult_hospital_beds_7_day_sum <dbl> ~
## $ all_adult_hospital_inpatient_beds_7_day_sum <dbl> ~
## $ inpatient_beds_used_7_day_sum <dbl> ~
## $ all_adult_hospital_inpatient_bed_occupied_7_day_sum <dbl> ~
## $ inpatient_beds_used_covid_7_day_sum <dbl> ~
## $ total_adult_patients_hospitalized_confirmed_and_suspected_covid_7_day_sum <dbl> ~
## $ total_adult_patients_hospitalized_confirmed_covid_7_day_sum <dbl> ~
## $ total_pediatric_patients_hospitalized_confirmed_and_suspected_covid_7_day_sum <dbl> ~
## $ total_pediatric_patients_hospitalized_confirmed_covid_7_day_sum <dbl> ~
## $ inpatient_beds_7_day_sum <dbl> ~
## $ total_icu_beds_7_day_sum <dbl> ~
## $ total_staffed_adult_icu_beds_7_day_sum <dbl> ~
## $ icu_beds_used_7_day_sum <dbl> ~
## $ staffed_adult_icu_bed_occupancy_7_day_sum <dbl> ~
## $ staffed_icu_adult_patients_confirmed_and_suspected_covid_7_day_sum <dbl> ~
## $ staffed_icu_adult_patients_confirmed_covid_7_day_sum <dbl> ~
## $ total_patients_hospitalized_confirmed_influenza_7_day_sum <dbl> ~
## $ icu_patients_confirmed_influenza_7_day_sum <dbl> ~
## $ total_patients_hospitalized_confirmed_influenza_and_covid_7_day_sum <dbl> ~
## $ total_beds_7_day_coverage <dbl> ~
## $ all_adult_hospital_beds_7_day_coverage <dbl> ~
## $ all_adult_hospital_inpatient_beds_7_day_coverage <dbl> ~
## $ inpatient_beds_used_7_day_coverage <dbl> ~
## $ all_adult_hospital_inpatient_bed_occupied_7_day_coverage <dbl> ~
## $ inpatient_beds_used_covid_7_day_coverage <dbl> ~
## $ total_adult_patients_hospitalized_confirmed_and_suspected_covid_7_day_coverage <dbl> ~
## $ total_adult_patients_hospitalized_confirmed_covid_7_day_coverage <dbl> ~
## $ total_pediatric_patients_hospitalized_confirmed_and_suspected_covid_7_day_coverage <dbl> ~
## $ total_pediatric_patients_hospitalized_confirmed_covid_7_day_coverage <dbl> ~
## $ inpatient_beds_7_day_coverage <dbl> ~
## $ total_icu_beds_7_day_coverage <dbl> ~
## $ total_staffed_adult_icu_beds_7_day_coverage <dbl> ~
## $ icu_beds_used_7_day_coverage <dbl> ~
## $ staffed_adult_icu_bed_occupancy_7_day_coverage <dbl> ~
## $ staffed_icu_adult_patients_confirmed_and_suspected_covid_7_day_coverage <dbl> ~
## $ staffed_icu_adult_patients_confirmed_covid_7_day_coverage <dbl> ~
## $ total_patients_hospitalized_confirmed_influenza_7_day_coverage <dbl> ~
## $ icu_patients_confirmed_influenza_7_day_coverage <dbl> ~
## $ total_patients_hospitalized_confirmed_influenza_and_covid_7_day_coverage <dbl> ~
## $ previous_day_admission_adult_covid_confirmed_7_day_sum <dbl> ~
## $ `previous_day_admission_adult_covid_confirmed_18-19_7_day_sum` <dbl> ~
## $ `previous_day_admission_adult_covid_confirmed_20-29_7_day_sum` <dbl> ~
## $ `previous_day_admission_adult_covid_confirmed_30-39_7_day_sum` <dbl> ~
## $ `previous_day_admission_adult_covid_confirmed_40-49_7_day_sum` <dbl> ~
## $ `previous_day_admission_adult_covid_confirmed_50-59_7_day_sum` <dbl> ~
## $ `previous_day_admission_adult_covid_confirmed_60-69_7_day_sum` <dbl> ~
## $ `previous_day_admission_adult_covid_confirmed_70-79_7_day_sum` <dbl> ~
## $ `previous_day_admission_adult_covid_confirmed_80+_7_day_sum` <dbl> ~
## $ previous_day_admission_adult_covid_confirmed_unknown_7_day_sum <dbl> ~
## $ previous_day_admission_pediatric_covid_confirmed_7_day_sum <dbl> ~
## $ previous_day_covid_ED_visits_7_day_sum <dbl> ~
## $ previous_day_admission_adult_covid_suspected_7_day_sum <dbl> ~
## $ `previous_day_admission_adult_covid_suspected_18-19_7_day_sum` <dbl> ~
## $ `previous_day_admission_adult_covid_suspected_20-29_7_day_sum` <dbl> ~
## $ `previous_day_admission_adult_covid_suspected_30-39_7_day_sum` <dbl> ~
## $ `previous_day_admission_adult_covid_suspected_40-49_7_day_sum` <dbl> ~
## $ `previous_day_admission_adult_covid_suspected_50-59_7_day_sum` <dbl> ~
## $ `previous_day_admission_adult_covid_suspected_60-69_7_day_sum` <dbl> ~
## $ `previous_day_admission_adult_covid_suspected_70-79_7_day_sum` <dbl> ~
## $ `previous_day_admission_adult_covid_suspected_80+_7_day_sum` <dbl> ~
## $ previous_day_admission_adult_covid_suspected_unknown_7_day_sum <dbl> ~
## $ previous_day_admission_pediatric_covid_suspected_7_day_sum <dbl> ~
## $ previous_day_total_ED_visits_7_day_sum <dbl> ~
## $ previous_day_admission_influenza_confirmed_7_day_sum <dbl> ~
## $ geocoded_hospital_address <chr> ~
## $ hhs_ids <chr> ~
## $ previous_day_admission_adult_covid_confirmed_7_day_coverage <dbl> ~
## $ previous_day_admission_pediatric_covid_confirmed_7_day_coverage <dbl> ~
## $ previous_day_admission_adult_covid_suspected_7_day_coverage <dbl> ~
## $ previous_day_admission_pediatric_covid_suspected_7_day_coverage <dbl> ~
## $ previous_week_personnel_covid_vaccinated_doses_administered_7_day <dbl> ~
## $ total_personnel_covid_vaccinated_doses_none_7_day <dbl> ~
## $ total_personnel_covid_vaccinated_doses_one_7_day <dbl> ~
## $ total_personnel_covid_vaccinated_doses_all_7_day <dbl> ~
## $ previous_week_patients_covid_vaccinated_doses_one_7_day <dbl> ~
## $ previous_week_patients_covid_vaccinated_doses_all_7_day <dbl> ~
## $ is_corrected <lgl> ~
##
## Hospital Subtype Counts:
## # A tibble: 4 x 2
## hospital_subtype n
## <chr> <int>
## 1 Childrens Hospitals 7503
## 2 Critical Access Hospitals 106952
## 3 Long Term 27474
## 4 Short Term 257934
##
## Records other than 50 states and DC
## # A tibble: 5 x 2
## state n
## <chr> <int>
## 1 AS 25
## 2 GU 160
## 3 MP 80
## 4 PR 4400
## 5 VI 160
##
## Record types for key metrics
## # A tibble: 8 x 5
## name `NA` Positive `Value -999999` Total
## <chr> <int> <int> <int> <int>
## 1 all_adult_hospital_beds_7_day_avg 11667 387469 727 399863
## 2 all_adult_hospital_inpatient_bed_occupi~ 3328 364400 32135 399863
## 3 icu_beds_used_7_day_avg 1649 350757 47457 399863
## 4 inpatient_beds_7_day_avg 1730 396567 1566 399863
## 5 staffed_icu_adult_patients_confirmed_an~ 4251 279744 115868 399863
## 6 total_adult_patients_hospitalized_confi~ 2372 278715 118776 399863
## 7 total_beds_7_day_avg 6632 392858 373 399863
## 8 total_icu_beds_7_day_avg 2064 377884 19915 399863## `stat_bin()` using `bins = 30`. Pick better value with `binwidth`.
saveToRDS(indivHosp_20220221, ovrWriteError=FALSE)The post-processing capabilities are included:
# Create pivoted burden data
burdenPivotList_220220 <- postProcessCDCDaily(cdc_daily_220220,
dataThruLabel="Jan 2022",
keyDatesBurden=c("2022-01-31", "2021-07-31",
"2021-01-31", "2020-07-31"
),
keyDatesVaccine=c("2021-12-31", "2021-09-30",
"2021-06-30", "2021-03-31"
),
returnData=TRUE
)## Joining, by = "state"##
## *** File has been checked for uniqueness by: state date name## Warning: Removed 24 row(s) containing missing values (geom_path).
## Warning: Removed 24 rows containing missing values (position_stack).
## Warning: Removed 24 rows containing missing values (position_stack).
## Warning: Removed 9 row(s) containing missing values (geom_path).
The hospital summaries are also added:
# Can be run only as-needed
dfStateAgeBucket2019 <- readPopStateAge("./RInputFiles/sc-est2019-agesex-civ.csv") %>%
filterPopStateAge(keyCol="POPEST2019_CIV", keyColName="pop2019") %>%
bucketPopStateAge(popVar="pop2019")##
## -- Column specification --------------------------------------------------------
## cols(
## SUMLEV = col_character(),
## REGION = col_double(),
## DIVISION = col_double(),
## STATE = col_double(),
## NAME = col_character(),
## SEX = col_double(),
## AGE = col_double(),
## ESTBASE2010_CIV = col_double(),
## POPEST2010_CIV = col_double(),
## POPEST2011_CIV = col_double(),
## POPEST2012_CIV = col_double(),
## POPEST2013_CIV = col_double(),
## POPEST2014_CIV = col_double(),
## POPEST2015_CIV = col_double(),
## POPEST2016_CIV = col_double(),
## POPEST2017_CIV = col_double(),
## POPEST2018_CIV = col_double(),
## POPEST2019_CIV = col_double()
## )##
## *** File has been checked for uniqueness by: NAME SEX AGE
##
## [1] TRUE
## [1] TRUE
## [1] TRUE
##
## PASSED CHECK: United States total is the sum of states and DC
##
##
## PASSED CHECK: Age 999 total is the sum of the ages
##
##
## PASSED CHECK: Sex 0 total is the sum of the sexes
# Create hospitalized per capita data
hospPerCap_220220 <- hospAgePerCapita(dfStateAgeBucket2019,
lst=burdenPivotList_220220,
popVar="pop2019",
excludeState=c(),
cumStartDate="2020-07-15"
)## Warning: Removed 18 row(s) containing missing values (geom_path).
The one-page CFR plot capability is included:
# Create CFR plots for select states
cfrStates <- list("FL"=list(keyState="FL", minDate="2020-08-01", multDeath=70),
"LA"=list(keyState="LA", minDate="2020-08-01", multDeath=80),
"CA"=list(keyState="CA", minDate="2020-08-01", multDeath=100),
"IL"=list(keyState="IL", minDate="2020-08-01", multDeath=100)
)
purrr::walk(cfrStates, .f=function(x) onePageCFRPlot(burdenPivotList_220220$dfPivot,
keyState=x$keyState,
minDate=x$minDate,
multDeath=x$multDeath
)
)
The peaks and valleys plots are included:
# Burden data
cdc_daily_220220$dfPerCapita %>%
mutate(regn=c(as.character(state.region), "South")[match(state, c(state.abb, "DC"))]) %>%
makePeakValley(numVar=c("new_deaths", "new_cases", "inp"),
windowWidth = 71,
rollMean=7,
facetVar=c("regn"),
fnNumVar=list("new_deaths"=function(x) x,
"new_cases"=function(x) x/1000,
"inp"=function(x) x/1000
),
fnPeak=list("new_deaths"=function(x) x+100,
"new_cases"=function(x) x+10,
"inp"=function(x) x+10
),
fnValley=list("new_deaths"=function(x) x-100,
"new_cases"=function(x) x-5,
"inp"=function(x) x-5
),
useTitle=c("new_deaths"="US coronavirus deaths",
"new_cases"="US coronavirus cases",
"inp"="US coronavirus total hospitalized"
),
yLab=c("new_deaths"="Rolling 7-day mean deaths",
"new_cases"="Rolling 7-day mean cases (000)",
"inp"="Rolling 7-day mean in hospital (000)"
)
)## Warning: Removed 6 row(s) containing missing values (geom_path).
## Warning: Removed 6 row(s) containing missing values (geom_path).
## Warning: Removed 6 row(s) containing missing values (geom_path).
## # A tibble: 3,107 x 11
## date regn new_deaths new_cases inp new_deaths_isPe~ new_cases_isPeak
## <date> <chr> <dbl> <dbl> <dbl> <lgl> <lgl>
## 1 2020-01-01 Nort~ NA NA NA FALSE FALSE
## 2 2020-01-01 South NA NA NA FALSE FALSE
## 3 2020-01-01 West NA NA NA FALSE FALSE
## 4 2020-01-02 Nort~ NA NA NA FALSE FALSE
## 5 2020-01-02 South NA NA NA FALSE FALSE
## 6 2020-01-02 West NA NA NA FALSE FALSE
## 7 2020-01-03 Nort~ NA NA NA FALSE FALSE
## 8 2020-01-03 South NA NA NA FALSE FALSE
## 9 2020-01-03 West NA NA NA FALSE FALSE
## 10 2020-01-04 Nort~ 0 0 0 FALSE FALSE
## # ... with 3,097 more rows, and 4 more variables: inp_isPeak <lgl>,
## # new_deaths_isValley <lgl>, new_cases_isValley <lgl>, inp_isValley <lgl># Vaccinations data for states with 8+ million population
cdc_daily_220220$dfPerCapita %>%
inner_join(getStateData(), by=c("state")) %>%
filter(pop >= 8000000) %>%
select(date, state, vxa, vxc) %>%
arrange(date, state) %>%
group_by(state) %>%
mutate(across(c(vxa, vxc), .fns=function(x) x-lag(x))) %>%
ungroup() %>%
mutate(regn=c(as.character(state.region), "South")[match(state, c(state.abb, "DC"))]) %>%
filter(date >= "2020-12-01") %>%
makePeakValley(numVar=c("vxc", "vxa"),
windowWidth = 29,
rollMean=21,
facetVar=c("state"),
fnNumVar=list("vxa"=function(x) x/1000,
"vxc"=function(x) x/1000
),
fnPeak=list("vxa"=function(x) x+25*max(x, na.rm=TRUE)/400,
"vxc"=function(x) x+25*max(x, na.rm=TRUE)/400
),
fnValley=list("vxa"=function(x) x-25*max(x, na.rm=TRUE)/400,
"vxc"=function(x) x-25*max(x, na.rm=TRUE)/400
),
fnGroupFacet=TRUE,
useTitle=c("vxa"="Vaccines adminsitered (US)",
"vxc"="Became fully vaccinated (US)"
),
yLab=c("vxa"="Rolling 21-day mean administered (000)",
"vxc"="Rolling 21-day mean completed (000)"
)
)## Warning: Removed 20 row(s) containing missing values (geom_path).
## Warning: Removed 20 row(s) containing missing values (geom_path).
## # A tibble: 5,364 x 8
## date state vxc vxa vxc_isPeak vxa_isPeak vxc_isValley vxa_isValley
## <date> <chr> <dbl> <dbl> <lgl> <lgl> <lgl> <lgl>
## 1 2020-12-01 CA NA NA FALSE FALSE FALSE FALSE
## 2 2020-12-01 FL NA NA FALSE FALSE FALSE FALSE
## 3 2020-12-01 GA NA NA FALSE FALSE FALSE FALSE
## 4 2020-12-01 IL NA NA FALSE FALSE FALSE FALSE
## 5 2020-12-01 MI NA NA FALSE FALSE FALSE FALSE
## 6 2020-12-01 NC NA NA FALSE FALSE FALSE FALSE
## 7 2020-12-01 NJ NA NA FALSE FALSE FALSE FALSE
## 8 2020-12-01 NY NA NA FALSE FALSE FALSE FALSE
## 9 2020-12-01 OH NA NA FALSE FALSE FALSE FALSE
## 10 2020-12-01 PA NA NA FALSE FALSE FALSE FALSE
## # ... with 5,354 more rowsThe hospital utlization plots are included:
indivHosp_20220221 %>%
filter(state %in% c(state.abb, "DC"),
collection_week==max(collection_week)
) %>%
pull(hospital_pk) %>%
plotHospitalUtilization(df=indivHosp_20220221, keyHosp=., plotTitle="US Hospitals Summed")
Imputed hospital utilization data are also created, using functional form:
# Impute values for hospital capacity
imputeNACapacity <- function(df,
keyStates=c(state.abb, "DC"),
varMapper=hhsMapper,
varsToImpute=c("total_beds", "adult_beds"),
varUsedToImpute=c("inpatient_beds")
) {
# FUNCTION ARGUMENTS:
# df: the initial data frame
# keyState: states to include for filtering
# varMapper: variables to include and output names (named vector of form c("original name"="modified name"))
# varsToImpute: variables to be imputed
# varUsedToImpute: percent changes in this variable assumed to drive percent changes in varsToImpute if NA
df %>%
filter(state %in% all_of(keyStates)) %>%
colSelector(c("state", "collection_week", "hospital_pk", names(varMapper))) %>%
colRenamer(varMapper) %>%
mutate(across(where(is.numeric), .fns=function(x) ifelse(is.na(x), NA, ifelse(x==-999999, NA, x)))) %>%
arrange(hospital_pk, collection_week) %>%
group_by(hospital_pk) %>%
mutate(across(all_of(varsToImpute),
.fns=function(x) testImputeNA(x=x, y=get(varUsedToImpute), naValues=-999999)
)
) %>%
group_by(state, collection_week) %>%
summarize(across(where(is.numeric), .fns=sum, na.rm=TRUE), n=n(),.groups="drop")
}
modStateHosp_20220221 <- imputeNACapacity(indivHosp_20220221)The function is split so that it is more generic:
# Select and filter as needed
skinnyHHS <- function(df,
keyStates=c(state.abb, "DC"),
idCols=c("state", "collection_week", "hospital_pk"),
varMapper=hhsMapper
) {
# FUNCTION ARGUMENTS:
# df: the initial data frame
# keyState: states to include for filtering
# varMapper: variables to include and output names (named vector of form c("original name"="modified name"))
df %>%
filter(state %in% all_of(keyStates)) %>%
colSelector(c(all_of(idCols), names(varMapper))) %>%
colRenamer(varMapper)
}
# Impute values for hospital capacity
imputeNACapacity <- function(df,
extraNA=c(-999999),
convertAllNA=TRUE,
idVars=c("hospital_pk"),
sortVars=c("collection_week"),
varsToImpute=c("total_beds", "adult_beds"),
varUsedToImpute=c("inpatient_beds")
) {
# FUNCTION ARGUMENTS:
# df: the initial data frame
# extraNA: values that should be treated as NA
# convertAllNA: boolean, should all extraNA values be converted in all numeric columns?
# if FALSE, extraNA values will not be converted, though imputing will treat as NA
# varsToImpute: variables to be imputed
# varUsedToImpute: percent changes in this variable assumed to drive percent changes in varsToImpute if NA
# Convert NA if requested
if(isTRUE(convertAllNA)) {
df <- df %>%
mutate(across(where(is.numeric),
.fns=function(x) ifelse(is.na(x), NA, ifelse(x %in% all_of(extraNA), NA, x))
)
)
}
# Impute values and return data
df %>%
arrange(across(all_of(c(idVars, sortVars)))) %>%
group_by(across(all_of(idVars))) %>%
mutate(across(all_of(varsToImpute),
.fns=function(x) testImputeNA(x=x, y=get(varUsedToImpute), naValues=extraNA)
)
) %>%
ungroup()
}
sumImputedHHS <- function(df,
groupVars=c("state", "collection_week")) {
# FUNCTION ARGUMENTS:
# df: the initial data frame
# groupVars: variables for summing the data to
df %>%
group_by(across(all_of(groupVars))) %>%
summarize(across(where(is.numeric), .fns=sum, na.rm=TRUE), n=n(),.groups="drop")
}
identical(skinnyHHS(indivHosp_20220221) %>%
imputeNACapacity() %>%
sumImputedHHS(),
modStateHosp_20220221
)## [1] TRUEUpdated maps with imputed capacity are created:
modStateHosp_20220221 <- skinnyHHS(indivHosp_20220221) %>%
imputeNACapacity() %>%
sumImputedHHS()
# ICU summary
createGeoMap(modStateHosp_20220221,
yVars=list("pctCovidICU"=c("label"="Covid", "color"="red"),
"pctICU"=c("label"="Total", "color"="black")
),
fullList=list("pctICU"=expression(icu_beds_occupied/icu_beds),
"pctCovidICU"=expression(adult_icu_covid/icu_beds)
),
plotTitle="Average % ICU Capacity Filled by Week",
plotSubtitle="August 2020 to January 2022",
plotScaleLabel="% ICU\nUsed",
returnData=FALSE
)
# Adult beds summary
createGeoMap(modStateHosp_20220221 %>% filter(!(state %in% c("CT", "DE", "SD", "AK"))),
yVars=list("pctCovidAdult"=c("label"="Covid", "color"="red"),
"pctAdult"=c("label"="Total", "color"="black")
),
fullList=list("pctAdult"=expression(adult_beds_occupied/adult_beds),
"pctCovidAdult"=expression(adult_beds_covid/adult_beds)
),
plotTitle="Average % Adult Beds Capacity Filled by Week",
plotSubtitle="August 2020 to January 2022\n(AK, CT, DE, and SD data excluded)",
plotScaleLabel="% Adult\nBeds\nUsed",
returnData=FALSE
)
Example Data Update
The function is run to download and process the latest data:
readList <- list("cdcDaily"="./RInputFiles/Coronavirus/CDC_dc_downloaded_220304.csv",
"cdcHosp"="./RInputFiles/Coronavirus/CDC_h_downloaded_220304.csv",
"vax"="./RInputFiles/Coronavirus/vaxData_downloaded_220304.csv"
)
compareList <- list("cdcDaily"=readFromRDS("cdc_daily_220220")$dfRaw$cdcDaily,
"cdcHosp"=readFromRDS("cdc_daily_220220")$dfRaw$cdcHosp,
"vax"=readFromRDS("cdc_daily_220220")$dfRaw$vax
)
cdc_daily_220304 <- readRunCDCDaily(thruLabel="Mar 2, 2022",
downloadTo=lapply(readList, FUN=function(x) if(file.exists(x)) NA else x),
readFrom=readList,
compareFile=compareList,
writeLog=NULL,
useClusters=readFromRDS("cdc_daily_210528")$useClusters,
weightedMeanAggs=c("tcpm7", "tdpm7", "cpm7", "dpm7", "hpm7",
"vxcpm7", "vxcgte65pct"
),
skipAssessmentPlots=FALSE,
brewPalette="Paired"
)##
## -- Column specification --------------------------------------------------------
## cols(
## submission_date = col_character(),
## state = col_character(),
## tot_cases = col_double(),
## conf_cases = col_double(),
## prob_cases = col_double(),
## new_case = col_double(),
## pnew_case = col_double(),
## tot_death = col_double(),
## conf_death = col_double(),
## prob_death = col_double(),
## new_death = col_double(),
## pnew_death = col_double(),
## created_at = col_character(),
## consent_cases = col_character(),
## consent_deaths = col_character()
## )##
## *** File has been checked for uniqueness by: state date
##
##
## Checking for similarity of: column names
## In reference but not in current:
## In current but not in reference:
##
## Checking for similarity of: date
## In reference but not in current: 0
## In current but not in reference: 12
##
## Checking for similarity of: state
## In reference but not in current:
## In current but not in reference:
##
##
## ***Differences of at least 5 and at least 5%
##
## date name newValue refValue absDelta pctDelta
## 1 2020-03-03 tot_cases 175 188 13 0.07162534
## 2 2022-02-13 new_deaths 615 446 169 0.31856739
## 3 2022-02-06 new_deaths 609 472 137 0.25346901
## 4 2022-02-12 new_deaths 891 695 196 0.24716267
## 5 2022-02-05 new_deaths 1158 989 169 0.15742897
## 6 2022-01-30 new_deaths 869 796 73 0.08768769
## 7 2022-02-07 new_deaths 3177 3000 177 0.05730937
## 8 2022-02-08 new_deaths 3704 3504 200 0.05549390
## 9 2022-02-03 new_deaths 2653 2515 138 0.05340557
## 10 2022-01-29 new_deaths 1469 1394 75 0.05239260
## 11 2022-02-11 new_deaths 2775 2638 137 0.05061888
## 12 2020-03-03 new_cases 51 64 13 0.22608696
## 13 2022-02-12 new_cases 66377 55089 11288 0.18586271
## 14 2022-02-13 new_cases 47803 40950 6853 0.15442858
## 15 2022-01-30 new_cases 155259 138089 17170 0.11706233
## 16 2022-02-05 new_cases 102256 91295 10961 0.11326214
## 17 2022-02-14 new_cases 178028 199342 21314 0.11296075
## 18 2022-02-11 new_cases 155537 172496 16959 0.10339813
## 19 2022-01-29 new_cases 215839 195076 20763 0.10105740
## 20 2020-03-07 new_cases 146 160 14 0.09150327
## 21 2021-10-31 new_cases 22766 20850 1916 0.08785767
## 22 2021-11-06 new_cases 32140 29452 2688 0.08728406
## 23 2021-10-24 new_cases 25952 23899 2053 0.08236545
## 24 2021-11-07 new_cases 28368 26379 1989 0.07266152
## 25 2020-03-06 new_cases 130 121 9 0.07171315
## 26 2021-10-23 new_cases 33628 31349 2279 0.07014790
## 27 2022-01-22 new_cases 320403 299989 20414 0.06581000
## 28 2022-02-06 new_cases 96184 90271 5913 0.06342549
## 29 2022-01-18 new_cases 861976 917498 55522 0.06240271
## 30 2020-03-09 new_cases 390 415 25 0.06211180
## 31 2022-01-31 new_cases 583405 620416 37011 0.06148921
## 32 2022-01-23 new_cases 310096 291779 18317 0.06086646
## 33 2021-11-14 new_cases 30649 28992 1657 0.05556580
## 34 2021-11-20 new_cases 42759 40531 2228 0.05349982
## 35 2021-12-25 new_cases 126095 119545 6550 0.05333008
## 36 2021-10-30 new_cases 31410 29822 1588 0.05186830
## 37 2021-05-24 new_cases 15400 16206 806 0.05100297
##
##
## ***Differences of at least 0 and at least 0.1%
##
## state name newValue refValue absDelta pctDelta
## 1 RI tot_deaths 1395734 1419362 23628 0.016786639
## 2 FL tot_deaths 22083833 22194156 110323 0.004983198
## 3 KY tot_deaths 4208985 4192009 16976 0.004041427
## 4 NC tot_deaths 7123584 7111646 11938 0.001677247
## 5 AL tot_deaths 6231481 6222060 9421 0.001512983
## 6 RI tot_cases 75530417 79533518 4003101 0.051631619
## 7 ME tot_cases 38751314 36951188 1800126 0.047557900
## 8 WA tot_cases 263650336 264272454 622118 0.002356852
## 9 KY tot_cases 270113183 269767213 345970 0.001281654
## 10 AL new_deaths 18381 17877 504 0.027800761
## 11 FL new_deaths 70406 68581 1825 0.026261449
## 12 WA new_deaths 11615 11316 299 0.026078235
## 13 KY new_deaths 13885 13565 320 0.023315118
## 14 NC new_deaths 22277 22148 129 0.005807541
## 15 ME new_cases 225203 212435 12768 0.058349595
## 16 RI new_cases 336543 354045 17502 0.050687240
## 17 WA new_cases 1396813 1410596 13783 0.009819018
## 18 KY new_cases 1265367 1258310 7057 0.005592633
## 19 SD new_cases 234285 234961 676 0.002881218
## 20 NC new_cases 2563976 2559793 4183 0.001632782
## 21 SC new_cases 1451483 1449247 2236 0.001541681
##
##
##
## Raw file for cdcDaily:
## Rows: 46,260
## Columns: 15
## $ date <date> 2021-03-11, 2021-02-12, 2021-03-01, 2020-02-04, 2020-0~
## $ state <chr> "KS", "UT", "CO", "AR", "AR", "CO", "PW", "UT", "MA", "~
## $ tot_cases <dbl> 297229, 359641, 438745, 0, 56199, 1222893, 0, 636992, 7~
## $ conf_cases <dbl> 241035, 359641, 411869, NA, NA, 1117524, NA, 636992, 65~
## $ prob_cases <dbl> 56194, 0, 26876, NA, NA, 105369, NA, 0, 45550, 321, NA,~
## $ new_cases <dbl> 0, 1060, 677, 0, 547, 6962, 0, 0, 451, 619, 69, 24010, ~
## $ pnew_case <dbl> 0, 0, 60, NA, 0, 1247, 0, 0, 46, 1, 10, 4196, 264, 3202~
## $ tot_deaths <dbl> 4851, 1785, 5952, 0, 674, 10953, 0, 3787, 17818, 805, 8~
## $ conf_death <dbl> NA, 1729, 5218, NA, NA, 9666, NA, 3635, 17458, 624, NA,~
## $ prob_death <dbl> NA, 56, 734, NA, NA, 1287, NA, 152, 360, 181, NA, NA, 1~
## $ new_deaths <dbl> 0, 11, 1, 0, 11, 20, 0, 0, 5, 3, 0, 345, 8, 190, 0, 3, ~
## $ pnew_death <dbl> 0, 2, 0, NA, 0, 4, 0, 0, 0, 0, 0, 0, 2, 0, 0, 0, 1, NA,~
## $ created_at <chr> "03/12/2021 03:20:13 PM", "02/13/2021 02:50:08 PM", "03~
## $ consent_cases <chr> "Agree", "Agree", "Agree", "Not agree", "Not agree", "A~
## $ consent_deaths <chr> "N/A", "Agree", "Agree", "Not agree", "Not agree", "Agr~##
## -- Column specification --------------------------------------------------------
## cols(
## .default = col_double(),
## state = col_character(),
## date = col_date(format = ""),
## geocoded_state = col_logical()
## )
## i Use `spec()` for the full column specifications.
##
## *** File has been checked for uniqueness by: state date
##
##
## Checking for similarity of: column names
## In reference but not in current:
## In current but not in reference:
##
## Checking for similarity of: date
## In reference but not in current: 0
## In current but not in reference: 11
##
## Checking for similarity of: state
## In reference but not in current:
## In current but not in reference:
##
##
## ***Differences of at least 5 and at least 5%
##
## date name newValue refValue absDelta pctDelta
## 1 2022-02-20 inp 58908 62620 3712 0.06108880
## 2 2022-02-20 hosp_adult 56750 60478 3728 0.06360255
##
##
## ***Differences of at least 0 and at least 0.1%
##
## state name newValue refValue absDelta pctDelta
## 1 ND inp 110569 109813 756 0.006860814
## 2 WV hosp_ped 4497 4703 206 0.044782609
## 3 ME hosp_ped 1594 1538 56 0.035759898
## 4 MA hosp_ped 10034 9724 310 0.031379694
## 5 IN hosp_ped 15429 15261 168 0.010948192
## 6 KY hosp_ped 15750 15913 163 0.010295929
## 7 VA hosp_ped 14705 14555 150 0.010252905
## 8 NJ hosp_ped 16251 16415 164 0.010041021
## 9 NV hosp_ped 4105 4067 38 0.009300049
## 10 SC hosp_ped 7732 7661 71 0.009224972
## 11 AL hosp_ped 17872 17976 104 0.005802276
## 12 VT hosp_ped 360 362 2 0.005540166
## 13 KS hosp_ped 4025 4005 20 0.004981320
## 14 NM hosp_ped 6428 6457 29 0.004501358
## 15 IA hosp_ped 6509 6481 28 0.004311008
## 16 NH hosp_ped 761 758 3 0.003949967
## 17 FL hosp_ped 82260 82509 249 0.003022413
## 18 TN hosp_ped 18633 18581 52 0.002794647
## 19 WY hosp_ped 784 786 2 0.002547771
## 20 CO hosp_ped 18126 18084 42 0.002319801
## 21 SD hosp_ped 3899 3891 8 0.002053915
## 22 GA hosp_ped 43658 43742 84 0.001922197
## 23 AR hosp_ped 10931 10911 20 0.001831334
## 24 UT hosp_ped 7634 7621 13 0.001704359
## 25 CT hosp_ped 5640 5649 9 0.001594472
## 26 HI hosp_ped 2016 2019 3 0.001486989
## 27 MS hosp_ped 9380 9368 12 0.001280137
## 28 AZ hosp_ped 23979 23949 30 0.001251878
## 29 IL hosp_ped 36164 36121 43 0.001189735
## 30 MN hosp_ped 13210 13224 14 0.001059242
## 31 ND hosp_adult 104829 102042 2787 0.026944328
##
##
##
## Raw file for cdcHosp:
## Rows: 39,269
## Columns: 117
## $ state <chr> ~
## $ date <date> ~
## $ critical_staffing_shortage_today_yes <dbl> ~
## $ critical_staffing_shortage_today_no <dbl> ~
## $ critical_staffing_shortage_today_not_reported <dbl> ~
## $ critical_staffing_shortage_anticipated_within_week_yes <dbl> ~
## $ critical_staffing_shortage_anticipated_within_week_no <dbl> ~
## $ critical_staffing_shortage_anticipated_within_week_not_reported <dbl> ~
## $ hospital_onset_covid <dbl> ~
## $ hospital_onset_covid_coverage <dbl> ~
## $ inpatient_beds <dbl> ~
## $ inpatient_beds_coverage <dbl> ~
## $ inpatient_beds_used <dbl> ~
## $ inpatient_beds_used_coverage <dbl> ~
## $ inp <dbl> ~
## $ inpatient_beds_used_covid_coverage <dbl> ~
## $ previous_day_admission_adult_covid_confirmed <dbl> ~
## $ previous_day_admission_adult_covid_confirmed_coverage <dbl> ~
## $ previous_day_admission_adult_covid_suspected <dbl> ~
## $ previous_day_admission_adult_covid_suspected_coverage <dbl> ~
## $ previous_day_admission_pediatric_covid_confirmed <dbl> ~
## $ previous_day_admission_pediatric_covid_confirmed_coverage <dbl> ~
## $ previous_day_admission_pediatric_covid_suspected <dbl> ~
## $ previous_day_admission_pediatric_covid_suspected_coverage <dbl> ~
## $ staffed_adult_icu_bed_occupancy <dbl> ~
## $ staffed_adult_icu_bed_occupancy_coverage <dbl> ~
## $ staffed_icu_adult_patients_confirmed_and_suspected_covid <dbl> ~
## $ staffed_icu_adult_patients_confirmed_and_suspected_covid_coverage <dbl> ~
## $ staffed_icu_adult_patients_confirmed_covid <dbl> ~
## $ staffed_icu_adult_patients_confirmed_covid_coverage <dbl> ~
## $ hosp_adult <dbl> ~
## $ total_adult_patients_hospitalized_confirmed_and_suspected_covid_coverage <dbl> ~
## $ total_adult_patients_hospitalized_confirmed_covid <dbl> ~
## $ total_adult_patients_hospitalized_confirmed_covid_coverage <dbl> ~
## $ hosp_ped <dbl> ~
## $ total_pediatric_patients_hospitalized_confirmed_and_suspected_covid_coverage <dbl> ~
## $ total_pediatric_patients_hospitalized_confirmed_covid <dbl> ~
## $ total_pediatric_patients_hospitalized_confirmed_covid_coverage <dbl> ~
## $ total_staffed_adult_icu_beds <dbl> ~
## $ total_staffed_adult_icu_beds_coverage <dbl> ~
## $ inpatient_beds_utilization <dbl> ~
## $ inpatient_beds_utilization_coverage <dbl> ~
## $ inpatient_beds_utilization_numerator <dbl> ~
## $ inpatient_beds_utilization_denominator <dbl> ~
## $ percent_of_inpatients_with_covid <dbl> ~
## $ percent_of_inpatients_with_covid_coverage <dbl> ~
## $ percent_of_inpatients_with_covid_numerator <dbl> ~
## $ percent_of_inpatients_with_covid_denominator <dbl> ~
## $ inpatient_bed_covid_utilization <dbl> ~
## $ inpatient_bed_covid_utilization_coverage <dbl> ~
## $ inpatient_bed_covid_utilization_numerator <dbl> ~
## $ inpatient_bed_covid_utilization_denominator <dbl> ~
## $ adult_icu_bed_covid_utilization <dbl> ~
## $ adult_icu_bed_covid_utilization_coverage <dbl> ~
## $ adult_icu_bed_covid_utilization_numerator <dbl> ~
## $ adult_icu_bed_covid_utilization_denominator <dbl> ~
## $ adult_icu_bed_utilization <dbl> ~
## $ adult_icu_bed_utilization_coverage <dbl> ~
## $ adult_icu_bed_utilization_numerator <dbl> ~
## $ adult_icu_bed_utilization_denominator <dbl> ~
## $ geocoded_state <lgl> ~
## $ `previous_day_admission_adult_covid_confirmed_18-19` <dbl> ~
## $ `previous_day_admission_adult_covid_confirmed_18-19_coverage` <dbl> ~
## $ `previous_day_admission_adult_covid_confirmed_20-29` <dbl> ~
## $ `previous_day_admission_adult_covid_confirmed_20-29_coverage` <dbl> ~
## $ `previous_day_admission_adult_covid_confirmed_30-39` <dbl> ~
## $ `previous_day_admission_adult_covid_confirmed_30-39_coverage` <dbl> ~
## $ `previous_day_admission_adult_covid_confirmed_40-49` <dbl> ~
## $ `previous_day_admission_adult_covid_confirmed_40-49_coverage` <dbl> ~
## $ `previous_day_admission_adult_covid_confirmed_50-59` <dbl> ~
## $ `previous_day_admission_adult_covid_confirmed_50-59_coverage` <dbl> ~
## $ `previous_day_admission_adult_covid_confirmed_60-69` <dbl> ~
## $ `previous_day_admission_adult_covid_confirmed_60-69_coverage` <dbl> ~
## $ `previous_day_admission_adult_covid_confirmed_70-79` <dbl> ~
## $ `previous_day_admission_adult_covid_confirmed_70-79_coverage` <dbl> ~
## $ `previous_day_admission_adult_covid_confirmed_80+` <dbl> ~
## $ `previous_day_admission_adult_covid_confirmed_80+_coverage` <dbl> ~
## $ previous_day_admission_adult_covid_confirmed_unknown <dbl> ~
## $ previous_day_admission_adult_covid_confirmed_unknown_coverage <dbl> ~
## $ `previous_day_admission_adult_covid_suspected_18-19` <dbl> ~
## $ `previous_day_admission_adult_covid_suspected_18-19_coverage` <dbl> ~
## $ `previous_day_admission_adult_covid_suspected_20-29` <dbl> ~
## $ `previous_day_admission_adult_covid_suspected_20-29_coverage` <dbl> ~
## $ `previous_day_admission_adult_covid_suspected_30-39` <dbl> ~
## $ `previous_day_admission_adult_covid_suspected_30-39_coverage` <dbl> ~
## $ `previous_day_admission_adult_covid_suspected_40-49` <dbl> ~
## $ `previous_day_admission_adult_covid_suspected_40-49_coverage` <dbl> ~
## $ `previous_day_admission_adult_covid_suspected_50-59` <dbl> ~
## $ `previous_day_admission_adult_covid_suspected_50-59_coverage` <dbl> ~
## $ `previous_day_admission_adult_covid_suspected_60-69` <dbl> ~
## $ `previous_day_admission_adult_covid_suspected_60-69_coverage` <dbl> ~
## $ `previous_day_admission_adult_covid_suspected_70-79` <dbl> ~
## $ `previous_day_admission_adult_covid_suspected_70-79_coverage` <dbl> ~
## $ `previous_day_admission_adult_covid_suspected_80+` <dbl> ~
## $ `previous_day_admission_adult_covid_suspected_80+_coverage` <dbl> ~
## $ previous_day_admission_adult_covid_suspected_unknown <dbl> ~
## $ previous_day_admission_adult_covid_suspected_unknown_coverage <dbl> ~
## $ deaths_covid <dbl> ~
## $ deaths_covid_coverage <dbl> ~
## $ on_hand_supply_therapeutic_a_casirivimab_imdevimab_courses <dbl> ~
## $ on_hand_supply_therapeutic_b_bamlanivimab_courses <dbl> ~
## $ on_hand_supply_therapeutic_c_bamlanivimab_etesevimab_courses <dbl> ~
## $ previous_week_therapeutic_a_casirivimab_imdevimab_courses_used <dbl> ~
## $ previous_week_therapeutic_b_bamlanivimab_courses_used <dbl> ~
## $ previous_week_therapeutic_c_bamlanivimab_etesevimab_courses_used <dbl> ~
## $ icu_patients_confirmed_influenza <dbl> ~
## $ icu_patients_confirmed_influenza_coverage <dbl> ~
## $ previous_day_admission_influenza_confirmed <dbl> ~
## $ previous_day_admission_influenza_confirmed_coverage <dbl> ~
## $ previous_day_deaths_covid_and_influenza <dbl> ~
## $ previous_day_deaths_covid_and_influenza_coverage <dbl> ~
## $ previous_day_deaths_influenza <dbl> ~
## $ previous_day_deaths_influenza_coverage <dbl> ~
## $ total_patients_hospitalized_confirmed_influenza <dbl> ~
## $ total_patients_hospitalized_confirmed_influenza_and_covid <dbl> ~
## $ total_patients_hospitalized_confirmed_influenza_and_covid_coverage <dbl> ~
## $ total_patients_hospitalized_confirmed_influenza_coverage <dbl> ~##
## -- Column specification --------------------------------------------------------
## cols(
## .default = col_double(),
## Date = col_character(),
## Location = col_character()
## )
## i Use `spec()` for the full column specifications.
##
## *** File has been checked for uniqueness by: state date
##
##
## Checking for similarity of: column names
## In reference but not in current:
## In current but not in reference:
##
## Checking for similarity of: date
## In reference but not in current: 0
## In current but not in reference: 12
##
## Checking for similarity of: state
## In reference but not in current:
## In current but not in reference:
##
##
## ***Differences of at least 1 and at least 1%
##
## [1] date name newValue refValue absDelta pctDelta
## <0 rows> (or 0-length row.names)##
##
## ***Differences of at least 0 and at least 0.1%
##
## [1] state name newValue refValue absDelta pctDelta
## <0 rows> (or 0-length row.names)
##
##
##
## Raw file for vax:
## Rows: 28,760
## Columns: 82
## $ date <date> 2022-03-03, 2022-03-03, 2022-0~
## $ MMWR_week <dbl> 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9~
## $ state <chr> "NE", "NC", "TX", "CA", "AL", "~
## $ Distributed <dbl> 3775510, 20928600, 58996495, 86~
## $ Distributed_Janssen <dbl> 149600, 917900, 2609300, 368570~
## $ Distributed_Moderna <dbl> 1331380, 7886660, 21192040, 307~
## $ Distributed_Pfizer <dbl> 2294530, 12124040, 35195155, 51~
## $ Distributed_Unk_Manuf <dbl> 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0~
## $ Dist_Per_100K <dbl> 195177, 199546, 203465, 217827,~
## $ Distributed_Per_100k_12Plus <dbl> 233430, 232915, 244787, 255803,~
## $ Distributed_Per_100k_18Plus <dbl> 258892, 255621, 273182, 281107,~
## $ Distributed_Per_100k_65Plus <dbl> 1208330, 1195170, 1579880, 1474~
## $ vxa <dbl> 3086667, 16146189, 44500682, 71~
## $ Administered_12Plus <dbl> 2984836, 15667506, 42929166, 69~
## $ Administered_18Plus <dbl> 2784330, 14709026, 39448235, 63~
## $ Administered_65Plus <dbl> 818837, 4254447, 8966187, 14693~
## $ Administered_Janssen <dbl> 93421, 510563, 1535569, 2278802~
## $ Administered_Moderna <dbl> 1110596, 6003546, 16331212, 267~
## $ Administered_Pfizer <dbl> 1876552, 9631145, 26629492, 426~
## $ Administered_Unk_Manuf <dbl> 6098, 935, 4409, 15243, 477, 21~
## $ Admin_Per_100k <dbl> 159566, 153948, 153472, 181390,~
## $ Admin_Per_100k_12Plus <dbl> 184544, 174364, 178121, 205442,~
## $ Admin_Per_100k_18Plus <dbl> 190925, 179655, 182664, 208987,~
## $ Admin_Per_100k_65Plus <dbl> 262063, 242959, 240108, 251683,~
## $ Recip_Administered <dbl> 3099534, 16045766, 43251399, 71~
## $ Administered_Dose1_Recip <dbl> 1343086, 8641769, 20646737, 323~
## $ Administered_Dose1_Pop_Pct <dbl> 69.4, 82.4, 71.2, 81.8, 61.9, 6~
## $ Administered_Dose1_Recip_12Plus <dbl> 1287672, 8369683, 19753177, 309~
## $ Administered_Dose1_Recip_12PlusPop_Pct <dbl> 79.6, 93.1, 82.0, 91.9, 71.1, 7~
## $ Administered_Dose1_Recip_18Plus <dbl> 1192672, 7856851, 17992902, 284~
## $ Administered_Dose1_Recip_18PlusPop_Pct <dbl> 81.8, 95.0, 83.3, 92.9, 73.8, 7~
## $ Administered_Dose1_Recip_65Plus <dbl> 306831, 2160730, 3621726, 59741~
## $ Administered_Dose1_Recip_65PlusPop_Pct <dbl> 95.0, 95.0, 95.0, 95.0, 95.0, 9~
## $ vxc <dbl> 1210303, 6231369, 17444705, 278~
## $ vxcpoppct <dbl> 62.6, 59.4, 60.2, 70.5, 50.3, 5~
## $ Series_Complete_12Plus <dbl> 1164264, 6032933, 16838920, 267~
## $ Series_Complete_12PlusPop_Pct <dbl> 72.0, 67.1, 69.9, 79.5, 58.0, 6~
## $ vxcgte18 <dbl> 1078877, 5641201, 15431387, 245~
## $ vxcgte18pct <dbl> 74.0, 68.9, 71.5, 80.3, 60.3, 6~
## $ vxcgte65 <dbl> 284820, 1501131, 3209108, 51872~
## $ vxcgte65pct <dbl> 91.2, 85.7, 85.9, 88.9, 81.2, 8~
## $ Series_Complete_Janssen <dbl> 87478, 478497, 1339798, 2070288~
## $ Series_Complete_Moderna <dbl> 412033, 2159277, 6030460, 95834~
## $ Series_Complete_Pfizer <dbl> 709218, 3593448, 10073548, 1620~
## $ Series_Complete_Unk_Manuf <dbl> 1574, 147, 899, 4865, 654, 375,~
## $ Series_Complete_Janssen_12Plus <dbl> 87453, 478469, 1339351, 2069676~
## $ Series_Complete_Moderna_12Plus <dbl> 411994, 2159160, 6029642, 95826~
## $ Series_Complete_Pfizer_12Plus <dbl> 663259, 3395158, 9469060, 15076~
## $ Series_Complete_Unk_Manuf_12Plus <dbl> 1558, 146, 867, 4805, 654, 373,~
## $ Series_Complete_Janssen_18Plus <dbl> 87386, 477036, 1337802, 2062395~
## $ Series_Complete_Moderna_18Plus <dbl> 411823, 2156126, 6025533, 95567~
## $ Series_Complete_Pfizer_18Plus <dbl> 578185, 3007897, 8067213, 12955~
## $ Series_Complete_Unk_Manuf_18Plus <dbl> 1483, 142, 839, 4500, 650, 366,~
## $ Series_Complete_Janssen_65Plus <dbl> 6942, 54470, 177453, 201180, 36~
## $ Series_Complete_Moderna_65Plus <dbl> 138470, 721440, 1524709, 262127~
## $ Series_Complete_Pfizer_65Plus <dbl> 138496, 725155, 1506622, 236331~
## $ Series_Complete_Unk_Manuf_65Plus <dbl> 912, 66, 324, 1458, 415, 213, 1~
## $ Additional_Doses <dbl> 585237, 1574890, 6257276, 13507~
## $ Additional_Doses_Vax_Pct <dbl> 48.4, 25.3, 35.9, 48.5, 34.4, 3~
## $ Additional_Doses_12Plus <dbl> 585134, 1574750, 6256853, 13506~
## $ Additional_Doses_12Plus_Vax_Pct <dbl> 50.3, 26.1, 37.2, 50.5, 34.9, 3~
## $ Additional_Doses_18Plus <dbl> 565479, 1527656, 6053696, 12967~
## $ Additional_Doses_18Plus_Vax_Pct <dbl> 52.4, 27.1, 39.2, 52.8, 36.3, 4~
## $ Additional_Doses_50Plus <dbl> 364760, 1031758, 3720211, 72429~
## $ Additional_Doses_50Plus_Vax_Pct <dbl> 65.5, 34.3, 52.0, 63.8, 47.2, 5~
## $ Additional_Doses_65Plus <dbl> 211563, 585640, 1945440, 368152~
## $ Additional_Doses_65Plus_Vax_Pct <dbl> 74.3, 39.0, 60.6, 71.0, 56.9, 6~
## $ Additional_Doses_Moderna <dbl> 229241, 693169, 2749043, 586417~
## $ Additional_Doses_Pfizer <dbl> 349233, 854191, 3412442, 742991~
## $ Additional_Doses_Janssen <dbl> 6431, 27508, 95589, 213372, 152~
## $ Additional_Doses_Unk_Manuf <dbl> 332, 22, 202, 522, 81, 490, 73,~
## $ Administered_Dose1_Recip_5Plus <dbl> 1342804, 8639422, 20641353, 323~
## $ Administered_Dose1_Recip_5PlusPop_Pct <dbl> 74.5, 87.5, 76.4, 87.0, 65.9, 7~
## $ Series_Complete_5Plus <dbl> 1210243, 6230533, 17443348, 278~
## $ Series_Complete_5PlusPop_Pct <dbl> 67.1, 63.1, 64.6, 75.0, 53.5, 5~
## $ Administered_5Plus <dbl> 3086314, 16143045, 44494008, 71~
## $ Admin_Per_100k_5Plus <dbl> 171126, 163419, 164762, 192973,~
## $ Distributed_Per_100k_5Plus <dbl> 209340, 211864, 218465, 231812,~
## $ Series_Complete_Moderna_5Plus <dbl> 412011, 2159234, 6029941, 95831~
## $ Series_Complete_Pfizer_5Plus <dbl> 709196, 3592673, 10072988, 1620~
## $ Series_Complete_Janssen_5Plus <dbl> 87463, 478480, 1339521, 2069915~
## $ Series_Complete_Unk_Manuf_5Plus <dbl> 1573, 146, 898, 4864, 654, 373,~
##
## Column sums before and after applying filtering rules:
## # A tibble: 3 x 6
## isType tot_cases tot_deaths new_cases new_deaths n
## <chr> <dbl> <dbl> <dbl> <dbl> <dbl>
## 1 before 1.93e+10 3.28e+8 7.85e+7 931901 45489
## 2 after 1.92e+10 3.27e+8 7.80e+7 927317 39321
## 3 pctchg 4.93e- 3 4.33e-3 6.98e-3 0.00492 0.136
##
##
## Processed for cdcDaily:
## Rows: 39,321
## Columns: 6
## $ date <date> 2021-03-11, 2021-02-12, 2021-03-01, 2020-02-04, 2020-08-22~
## $ state <chr> "KS", "UT", "CO", "AR", "AR", "CO", "UT", "MA", "HI", "TX",~
## $ tot_cases <dbl> 297229, 359641, 438745, 0, 56199, 1222893, 636992, 704796, ~
## $ tot_deaths <dbl> 4851, 1785, 5952, 0, 674, 10953, 3787, 17818, 883, 33124, 7~
## $ new_cases <dbl> 0, 1060, 677, 0, 547, 6962, 0, 451, 69, 24010, 1028, 18811,~
## $ new_deaths <dbl> 0, 11, 1, 0, 11, 20, 0, 5, 0, 345, 8, 190, 3, 15, 7, 8, 0, ~
##
## Column sums before and after applying filtering rules:
## # A tibble: 3 x 5
## isType inp hosp_adult hosp_ped n
## <chr> <dbl> <dbl> <dbl> <dbl>
## 1 before 4.57e+7 3.93e+7 965351 39269
## 2 after 4.54e+7 3.91e+7 947960 37644
## 3 pctchg 4.81e-3 4.60e-3 0.0180 0.0414
##
##
## Processed for cdcHosp:
## Rows: 37,644
## Columns: 5
## $ date <date> 2020-10-18, 2020-10-13, 2020-10-12, 2020-10-08, 2020-10-06~
## $ state <chr> "VT", "NH", "ID", "MT", "HI", "NH", "NC", "DC", "MA", "MT",~
## $ inp <dbl> 2, 34, 221, 262, 124, 48, 1283, 156, 354, 207, 116, 102, 39~
## $ hosp_adult <dbl> 2, 34, 219, 259, 124, 48, 1246, 141, 347, 206, 109, 101, 38~
## $ hosp_ped <dbl> 0, 0, 2, 3, 0, 0, 34, 15, 7, 1, 3, 1, 10, 0, 0, 1, 6, 6, 7,~
##
## Column sums before and after applying filtering rules:
## # A tibble: 3 x 9
## isType vxa vxc vxcpoppct vxcgte65 vxcgte65pct vxcgte18 vxcgte18pct
## <chr> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl>
## 1 before 2.79e+11 1.18e+11 1050491. 3.15e+10 1622354. 1.11e+11 1255957.
## 2 after 1.34e+11 5.72e+10 882060. 1.52e+10 1450421 5.37e+10 1065505.
## 3 pctchg 5.20e- 1 5.16e- 1 0.160 5.16e- 1 0.106 5.17e- 1 0.152
## # ... with 1 more variable: n <dbl>
##
##
## Processed for vax:
## Rows: 22,695
## Columns: 9
## $ date <date> 2022-03-03, 2022-03-03, 2022-03-03, 2022-03-03, 2022-03-0~
## $ state <chr> "NE", "NC", "TX", "CA", "AL", "SC", "WV", "MN", "CO", "KS"~
## $ vxa <dbl> 3086667, 16146189, 44500682, 71671126, 6108052, 7287794, 2~
## $ vxc <dbl> 1210303, 6231369, 17444705, 27867605, 2466221, 2880832, 10~
## $ vxcpoppct <dbl> 62.6, 59.4, 60.2, 70.5, 50.3, 56.0, 56.8, 68.3, 69.3, 60.3~
## $ vxcgte65 <dbl> 284820, 1501131, 3209108, 5187220, 689667, 807207, 306687,~
## $ vxcgte65pct <dbl> 91.2, 85.7, 85.9, 88.9, 81.2, 86.1, 83.6, 95.0, 92.0, 89.5~
## $ vxcgte18 <dbl> 1078877, 5641201, 15431387, 24579521, 2300814, 2646739, 94~
## $ vxcgte18pct <dbl> 74.0, 68.9, 71.5, 80.3, 60.3, 65.6, 66.0, 78.2, 79.1, 71.1~
##
## Integrated per capita data file:
## Rows: 39,534
## Columns: 34
## $ date <date> 2020-01-01, 2020-01-01, 2020-01-01, 2020-01-01, 2020-01-0~
## $ state <chr> "AL", "HI", "IN", "LA", "MN", "MT", "NC", "TX", "AL", "HI"~
## $ tot_cases <dbl> NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA~
## $ tot_deaths <dbl> NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA~
## $ new_cases <dbl> NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA~
## $ new_deaths <dbl> NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA~
## $ inp <dbl> NA, 0, 0, NA, 0, 0, 0, 0, NA, 0, 0, NA, 0, 0, 0, 1877, 0, ~
## $ hosp_adult <dbl> NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA~
## $ hosp_ped <dbl> NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA~
## $ vxa <dbl> NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA~
## $ vxc <dbl> NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA~
## $ vxcpoppct <dbl> NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA~
## $ vxcgte65 <dbl> NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA~
## $ vxcgte65pct <dbl> NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA~
## $ vxcgte18 <dbl> NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA~
## $ vxcgte18pct <dbl> NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA~
## $ tcpm <dbl> NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA~
## $ tdpm <dbl> NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA~
## $ cpm <dbl> NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA~
## $ dpm <dbl> NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA~
## $ hpm <dbl> NA, 0.0000, 0.0000, NA, 0.0000, 0.0000, 0.0000, 0.0000, NA~
## $ ahpm <dbl> NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA~
## $ phpm <dbl> NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA~
## $ vxapm <dbl> NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA~
## $ vxcpm <dbl> NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA~
## $ tcpm7 <dbl> NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA~
## $ tdpm7 <dbl> NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA~
## $ cpm7 <dbl> NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA~
## $ dpm7 <dbl> NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA~
## $ hpm7 <dbl> NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA~
## $ ahpm7 <dbl> NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA~
## $ phpm7 <dbl> NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA~
## $ vxapm7 <dbl> NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA~
## $ vxcpm7 <dbl> NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA~## Warning in showSRID(uprojargs, format = "PROJ", multiline = "NO", prefer_proj =
## prefer_proj): Discarded datum unknown in CRS definition
saveToRDS(cdc_daily_220304, ovrWriteError=FALSE)
# Run for latest data, save as RDS
indivHosp_20220304 <- downloadReadHospitalData(loc="./RInputFiles/Coronavirus/HHS_Hospital_20220304.csv")##
## -- Column specification --------------------------------------------------------
## cols(
## .default = col_double(),
## hospital_pk = col_character(),
## collection_week = col_date(format = ""),
## state = col_character(),
## ccn = col_character(),
## hospital_name = col_character(),
## address = col_character(),
## city = col_character(),
## zip = col_character(),
## hospital_subtype = col_character(),
## fips_code = col_character(),
## is_metro_micro = col_logical(),
## geocoded_hospital_address = col_character(),
## hhs_ids = col_character(),
## is_corrected = col_logical()
## )
## i Use `spec()` for the full column specifications.## Rows: 409,797
## Columns: 109
## $ hospital_pk <chr> ~
## $ collection_week <date> ~
## $ state <chr> ~
## $ ccn <chr> ~
## $ hospital_name <chr> ~
## $ address <chr> ~
## $ city <chr> ~
## $ zip <chr> ~
## $ hospital_subtype <chr> ~
## $ fips_code <chr> ~
## $ is_metro_micro <lgl> ~
## $ total_beds_7_day_avg <dbl> ~
## $ all_adult_hospital_beds_7_day_avg <dbl> ~
## $ all_adult_hospital_inpatient_beds_7_day_avg <dbl> ~
## $ inpatient_beds_used_7_day_avg <dbl> ~
## $ all_adult_hospital_inpatient_bed_occupied_7_day_avg <dbl> ~
## $ inpatient_beds_used_covid_7_day_avg <dbl> ~
## $ total_adult_patients_hospitalized_confirmed_and_suspected_covid_7_day_avg <dbl> ~
## $ total_adult_patients_hospitalized_confirmed_covid_7_day_avg <dbl> ~
## $ total_pediatric_patients_hospitalized_confirmed_and_suspected_covid_7_day_avg <dbl> ~
## $ total_pediatric_patients_hospitalized_confirmed_covid_7_day_avg <dbl> ~
## $ inpatient_beds_7_day_avg <dbl> ~
## $ total_icu_beds_7_day_avg <dbl> ~
## $ total_staffed_adult_icu_beds_7_day_avg <dbl> ~
## $ icu_beds_used_7_day_avg <dbl> ~
## $ staffed_adult_icu_bed_occupancy_7_day_avg <dbl> ~
## $ staffed_icu_adult_patients_confirmed_and_suspected_covid_7_day_avg <dbl> ~
## $ staffed_icu_adult_patients_confirmed_covid_7_day_avg <dbl> ~
## $ total_patients_hospitalized_confirmed_influenza_7_day_avg <dbl> ~
## $ icu_patients_confirmed_influenza_7_day_avg <dbl> ~
## $ total_patients_hospitalized_confirmed_influenza_and_covid_7_day_avg <dbl> ~
## $ total_beds_7_day_sum <dbl> ~
## $ all_adult_hospital_beds_7_day_sum <dbl> ~
## $ all_adult_hospital_inpatient_beds_7_day_sum <dbl> ~
## $ inpatient_beds_used_7_day_sum <dbl> ~
## $ all_adult_hospital_inpatient_bed_occupied_7_day_sum <dbl> ~
## $ inpatient_beds_used_covid_7_day_sum <dbl> ~
## $ total_adult_patients_hospitalized_confirmed_and_suspected_covid_7_day_sum <dbl> ~
## $ total_adult_patients_hospitalized_confirmed_covid_7_day_sum <dbl> ~
## $ total_pediatric_patients_hospitalized_confirmed_and_suspected_covid_7_day_sum <dbl> ~
## $ total_pediatric_patients_hospitalized_confirmed_covid_7_day_sum <dbl> ~
## $ inpatient_beds_7_day_sum <dbl> ~
## $ total_icu_beds_7_day_sum <dbl> ~
## $ total_staffed_adult_icu_beds_7_day_sum <dbl> ~
## $ icu_beds_used_7_day_sum <dbl> ~
## $ staffed_adult_icu_bed_occupancy_7_day_sum <dbl> ~
## $ staffed_icu_adult_patients_confirmed_and_suspected_covid_7_day_sum <dbl> ~
## $ staffed_icu_adult_patients_confirmed_covid_7_day_sum <dbl> ~
## $ total_patients_hospitalized_confirmed_influenza_7_day_sum <dbl> ~
## $ icu_patients_confirmed_influenza_7_day_sum <dbl> ~
## $ total_patients_hospitalized_confirmed_influenza_and_covid_7_day_sum <dbl> ~
## $ total_beds_7_day_coverage <dbl> ~
## $ all_adult_hospital_beds_7_day_coverage <dbl> ~
## $ all_adult_hospital_inpatient_beds_7_day_coverage <dbl> ~
## $ inpatient_beds_used_7_day_coverage <dbl> ~
## $ all_adult_hospital_inpatient_bed_occupied_7_day_coverage <dbl> ~
## $ inpatient_beds_used_covid_7_day_coverage <dbl> ~
## $ total_adult_patients_hospitalized_confirmed_and_suspected_covid_7_day_coverage <dbl> ~
## $ total_adult_patients_hospitalized_confirmed_covid_7_day_coverage <dbl> ~
## $ total_pediatric_patients_hospitalized_confirmed_and_suspected_covid_7_day_coverage <dbl> ~
## $ total_pediatric_patients_hospitalized_confirmed_covid_7_day_coverage <dbl> ~
## $ inpatient_beds_7_day_coverage <dbl> ~
## $ total_icu_beds_7_day_coverage <dbl> ~
## $ total_staffed_adult_icu_beds_7_day_coverage <dbl> ~
## $ icu_beds_used_7_day_coverage <dbl> ~
## $ staffed_adult_icu_bed_occupancy_7_day_coverage <dbl> ~
## $ staffed_icu_adult_patients_confirmed_and_suspected_covid_7_day_coverage <dbl> ~
## $ staffed_icu_adult_patients_confirmed_covid_7_day_coverage <dbl> ~
## $ total_patients_hospitalized_confirmed_influenza_7_day_coverage <dbl> ~
## $ icu_patients_confirmed_influenza_7_day_coverage <dbl> ~
## $ total_patients_hospitalized_confirmed_influenza_and_covid_7_day_coverage <dbl> ~
## $ previous_day_admission_adult_covid_confirmed_7_day_sum <dbl> ~
## $ `previous_day_admission_adult_covid_confirmed_18-19_7_day_sum` <dbl> ~
## $ `previous_day_admission_adult_covid_confirmed_20-29_7_day_sum` <dbl> ~
## $ `previous_day_admission_adult_covid_confirmed_30-39_7_day_sum` <dbl> ~
## $ `previous_day_admission_adult_covid_confirmed_40-49_7_day_sum` <dbl> ~
## $ `previous_day_admission_adult_covid_confirmed_50-59_7_day_sum` <dbl> ~
## $ `previous_day_admission_adult_covid_confirmed_60-69_7_day_sum` <dbl> ~
## $ `previous_day_admission_adult_covid_confirmed_70-79_7_day_sum` <dbl> ~
## $ `previous_day_admission_adult_covid_confirmed_80+_7_day_sum` <dbl> ~
## $ previous_day_admission_adult_covid_confirmed_unknown_7_day_sum <dbl> ~
## $ previous_day_admission_pediatric_covid_confirmed_7_day_sum <dbl> ~
## $ previous_day_covid_ED_visits_7_day_sum <dbl> ~
## $ previous_day_admission_adult_covid_suspected_7_day_sum <dbl> ~
## $ `previous_day_admission_adult_covid_suspected_18-19_7_day_sum` <dbl> ~
## $ `previous_day_admission_adult_covid_suspected_20-29_7_day_sum` <dbl> ~
## $ `previous_day_admission_adult_covid_suspected_30-39_7_day_sum` <dbl> ~
## $ `previous_day_admission_adult_covid_suspected_40-49_7_day_sum` <dbl> ~
## $ `previous_day_admission_adult_covid_suspected_50-59_7_day_sum` <dbl> ~
## $ `previous_day_admission_adult_covid_suspected_60-69_7_day_sum` <dbl> ~
## $ `previous_day_admission_adult_covid_suspected_70-79_7_day_sum` <dbl> ~
## $ `previous_day_admission_adult_covid_suspected_80+_7_day_sum` <dbl> ~
## $ previous_day_admission_adult_covid_suspected_unknown_7_day_sum <dbl> ~
## $ previous_day_admission_pediatric_covid_suspected_7_day_sum <dbl> ~
## $ previous_day_total_ED_visits_7_day_sum <dbl> ~
## $ previous_day_admission_influenza_confirmed_7_day_sum <dbl> ~
## $ geocoded_hospital_address <chr> ~
## $ hhs_ids <chr> ~
## $ previous_day_admission_adult_covid_confirmed_7_day_coverage <dbl> ~
## $ previous_day_admission_pediatric_covid_confirmed_7_day_coverage <dbl> ~
## $ previous_day_admission_adult_covid_suspected_7_day_coverage <dbl> ~
## $ previous_day_admission_pediatric_covid_suspected_7_day_coverage <dbl> ~
## $ previous_week_personnel_covid_vaccinated_doses_administered_7_day <dbl> ~
## $ total_personnel_covid_vaccinated_doses_none_7_day <dbl> ~
## $ total_personnel_covid_vaccinated_doses_one_7_day <dbl> ~
## $ total_personnel_covid_vaccinated_doses_all_7_day <dbl> ~
## $ previous_week_patients_covid_vaccinated_doses_one_7_day <dbl> ~
## $ previous_week_patients_covid_vaccinated_doses_all_7_day <dbl> ~
## $ is_corrected <lgl> ~
##
## Hospital Subtype Counts:
## # A tibble: 4 x 2
## hospital_subtype n
## <chr> <int>
## 1 Childrens Hospitals 7690
## 2 Critical Access Hospitals 109641
## 3 Long Term 28161
## 4 Short Term 264305
##
## Records other than 50 states and DC
## # A tibble: 5 x 2
## state n
## <chr> <int>
## 1 AS 27
## 2 GU 164
## 3 MP 82
## 4 PR 4506
## 5 VI 164
##
## Record types for key metrics
## # A tibble: 8 x 5
## name `NA` Positive `Value -999999` Total
## <chr> <int> <int> <int> <int>
## 1 all_adult_hospital_beds_7_day_avg 15604 393445 748 409797
## 2 all_adult_hospital_inpatient_bed_occupi~ 3318 373556 32923 409797
## 3 icu_beds_used_7_day_avg 1649 359635 48513 409797
## 4 inpatient_beds_7_day_avg 1730 406462 1605 409797
## 5 staffed_icu_adult_patients_confirmed_an~ 4241 286438 119118 409797
## 6 total_adult_patients_hospitalized_confi~ 2362 285557 121878 409797
## 7 total_beds_7_day_avg 10392 399022 383 409797
## 8 total_icu_beds_7_day_avg 2064 387368 20365 409797## `stat_bin()` using `bins = 30`. Pick better value with `binwidth`.
saveToRDS(indivHosp_20220304, ovrWriteError=FALSE)Post-processing is run, including hospital summaries:
# Create pivoted burden data
burdenPivotList_220304 <- postProcessCDCDaily(cdc_daily_220304,
dataThruLabel="Feb 2022",
keyDatesBurden=c("2022-02-28", "2021-08-31",
"2021-02-28", "2020-08-31"
),
keyDatesVaccine=c("2022-02-28", "2021-10-31",
"2021-06-30", "2021-02-28"
),
returnData=TRUE
)## Joining, by = "state"##
## *** File has been checked for uniqueness by: state date name## Warning: Removed 24 row(s) containing missing values (geom_path).
## Warning: Removed 24 rows containing missing values (position_stack).
## Warning: Removed 24 rows containing missing values (position_stack).
## Warning: Removed 9 row(s) containing missing values (geom_path).
# Create hospitalized per capita data
hospPerCap_220304 <- hospAgePerCapita(readFromRDS("dfStateAgeBucket2019"),
lst=burdenPivotList_220304,
popVar="pop2019",
excludeState=c(),
cumStartDate="2020-07-15"
)## Warning: Removed 18 row(s) containing missing values (geom_path).
